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1. Introduce ion and Suimnary 

The purpose of this study was to evaluate. In detail the performance capabilities 
of specific SHUTTLE-based laser ranging systems, determine interface and support 
requirements, and to generate a preliminary design of a SHUTTLE-borne laser 
ranging experiment, Tlie study was conducted in two phases, referred to as 
(1) experiment definition and (2) experiment design. 

The goal of the experiment definition phase was to select, from the various 
laser ranging concepts, a viable approach to be used in the subsequent experi- 
ment design phase. The experiment objective Is to make laser ranging measure- 
ments to targets on the ground, from a Shuttle based laser ranging system. 

These measurements are to be used to determine the relative positions of the 
targets with respect to each other with an overall three dimensional recon- 
struction accuracy hatter than 2 cm ims. 'two basic approaches were identified. 
The first approach was to use relatively broad beamwidths, such thac each 
emitted laser pulse Illuminated all of the targets In a subset of the target 
grouping, referred to as a target grid. Thus, each pulse period would result 
in nearly simultaneous measurements to 2 or more targets. This approach mini- 
mises, but does not evade, dependence on orbital mechanics to reduce the data. 

It was found, however, that this approach was simply not a viable concept; the 
energy management problem was excessively severe for target grid dlmens ions of 
the order of magnitude needed to satisfy typical measurement objectives. 

The selected approach is to nuike the ranging measurements sequentially. However, 
in order to reduce the data. It is necessary to employ orbital meclumics tech- 
niques to Interpolate. It was unknown, at the start of this study, just how 
tightly the measurements had to be interlaced, to meet the grid reconstmctlon 
accuracy goal. Consequently, a mission simulation was developed to explore the 
effects of various measurement strategies on the experiment design. 

The stmulatioh csnployed a co-variance mfitrix analysis technique to determine 
the grid reconstruction accuracy, A Kalman filter was employed to incorporate 
tlie meashr omenta . Several measurement strategies were evaluated and a simple; 
strategy was selected, : A number of factors were found to be significant in 
accomplishing the desired reconstruction accuracy goals. First, we found that 
a Shuttle ephemeris errbr term coupled Intd the eo-varlance matrix of the target 
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grid position location unless widely spaced out-rigger targets were Included 
in the experiments Second^ vs found that at least two measurement passes were 
needed to reduce the reconstruction errors In all three axes. The second pass 
should ho as nearly orchagonal to the first pass as possible. We also explored 
mission phasing and found that unless the pass was within a few degrees of 
longitude of passing directly over the target grtd^ the recotistruction accuracy 
degraded considerably. Finally, we found that the interlace requirements were 
not severe. A measurement strategy which resulted in a grid survey in the 
approach period, a grid survey near the point of closest approach, and a final 
grid survey in the latter portion of the viewing opportunity provided excellent 
results. 

The modest interlace requirements do not require a high level of angular agility 
to satisfy the slewing and tracking requirements for the 9-target grid eval- 
uated in the mission simulation. Thus, it was postulated that one of the stable 
platforms which are candidates for the role of Shuttle/Spacelab general purpose 
experiment pointii\g systems could be employed to point the laser ranging experi- 
ment transmit and receive telescopes toward the targets in the grid. Tlie Small 
Instrument Pointing System (SIPS) was considered to be best suited for this 
purpose. Evaluation of the capabilities of the SIPS for the laser ranging exper- 
iment did not result in a clear decision between the SIPS mounted experiment and 
an experiment configured with a gimballed beam steering mirror. 


The study was therefore extended to include both coiicepts in the preliminary 
experiment deslgtx phase. The laser ranging systems for the two experiments were 
selected to be as much alike as possible. The results of the prellminavy design 
activities are discussed in Section 3. Both designs are considered implementable 
with state-of-tlic-art components and fully capable of meeting the experiment 
Ob.iectives, 
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2. Experiment Definition 

The first phase of the study was devoted to deterntlning constraints and 
requirements needed to define a viable experiment. A number of concepts and 
alternatives were explored and are reported herein. 

The experiment objective is to employ laser range measurements, made from the 
Shuttle borne laser ranging system to targets on the ground, to determine the 
relative spacings between these targets. The targets were assumed to be located 
roughly in a grid structure, with target arrays located at the corners of a 
square with one target In the middle. The grid dimension was assumed to be on 
the order of 10 km on a side, although the effect of varying grid size was 
considered. Each target in the grid was assumed to be composed of cube- 
corner retro-ref lectors, arranged such that the effective target radar cross- 

7 

section was on the order of 10 square meters. Grids of this sort could be used 
in a number of applications for various scientific or engineering purposes. 

One of the most fundamental questions to be addressed was the mensurement 
strategy to be employed. Two concepts were considered. The first concept was 
to employ s 1^1 taneous range n^asurements to two or more targets to minimize 
the dependency of the experiment on the Shuttle ephewerise The second concept was to 
employ sequential measurements to the targets, and to reduce the data with more 
sophisticated signal processing. The first concept is particularly appealing 
since differential range measurements could be used for data processing, 
eliminating a number of possible bias error sources. However, as discussed in 
Section 2.1, simultaneous range measuiements present a virtually Insurmountable 
energy management problem for modest to large grid dlmensidns. The second 
alternative, sequential range measurements requires sophisticated signal 
processing to extract the desired target relative location data. The range 
measurements need to be interlaced, i.e., measurements to each target in the 
grid must be made a number of times at various look angles. One possible 
measurement policy is to move from target to target during each interpulse 
period, resulting In a maximally interlaced data set. The problem with this 
measurement policy is that excessive slew rates are required to accomplish this 
for a modest grid size. The alternative is to dwell oii each target for a short 
period of time prior to moving to the next target. This policy mtiximlzes the 
data obtained when finite slew times and settling times are considered, bnt 
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the medsurementd are not maximally Interlaced. The effect of this measurement 
policy, or any other for that matter, on the achievable relative location 
accuracy can only be determined by simulation. This subject Is addressed in 
Section 2.2; a simple measurement policy was slwvn\ to be quite effective, and 
to not require excessive angular agility to accomplish the experiment objec- 
tives . 

A number of factors contribute to errors in pointing the transmitter and 
receiver at a specific target. Tltese fall Into two categories, location errors 
and attitude reference errors. Location errors include uncertainty as to the 
target location and uncertainties as to where the Shuttle really is. 

The pointing error attributable to these sources can be reduced to negligible 
magnitude once several ranging measurements are used to update the experiment 
state vector. Pointing errors attributable to attitude reference errors tend 
to be neatly stationary if the Shuttle is maintained in an inertially fixed 
attitude during the minutes of a tanging pass. Consequently, the major 
pointing problem is simply to find the first target, Once the target has been 
detected and a few range measurements are made, the system can complete the 
ranging pass with open- loop pointing. Section 2,3 describes our analyses of 
the acquisition process devised to cope with the initial pointing uncertaiiities , 

A liumbcr of stable platforms have been postulated to support various sol.ar^ 
and stellar astronomy instruments for Shuttle/Spacelab experiments. The con- 
cepts differ considerably, but are inte;\ded to permit experimenters to accom- 
plish their experiments using a general purpose, high stability, pointing system. 
Section 2,4 describes four of the concepts livvestlgated, and shows that the 
Small Instrument Pointing System (SIPS) appears to be the most suitable geiteral 
purpose poiitting system for the laser ranging experiment. 




2,1 Link Analysis 

The performance of the laser ranging experiment is dependent on the laser 
energy, the pulsewidth, the target radar cross-section, t}\e transmit beam- 


width, the receiver aperture, the link geometry, the various losses encountered, 
the signal processing techniques employed, and the background level. In 
general, it is possible to select the besnwidth, the transmit energy, the 
receiver aperture, and the target radar cross-taectioi^ to obtain suf ficient 
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return signal energy for confident range measurements for any particular link. 
There are a few limits, however, which affect these selections and the link 
feasibility. It is desirable to limit the peak transmitted energy density on 
the ground to values which are well below the accepted safety standards. Thus, 
for any given laser energy level, there Is a minimum transmit beamwidth. 

The optical receiver must employ relatively small active area detectors to 
obtain sufficient bandwidth for ranging purposes. Thus, there is a finite 
limit on the maximum receiver aperture which can be employed for any particular 
receiver fleld-of-view (FOV) . 

The first step in the link analysis process is to quantify these limits, and 
define a maximum performance system. Comparison of the achieved performance 
with the experiment objectives will then yield a performance margin, i.e. , the 
amount by Which the performance can be degraded and still achieve the experiment 
objectives. The final step in this process is to allocate the performance 
margin to system parameters in a maximally effective manner. 

2.1.1 Maximum Performance System Analysis 

The ‘first limit to be considered is the minimum transmit beantwidth. This limit 
is Imposed to assure rhat the maximum possible transmit energy density at the 
ground (E^) does not exceed the safe level for human exposure. Neglecting 
system losses, and assuming a Gaussian intensity profile, it can be shown that. 



■ ,2 

where “ full planar beamwidth (e power points) 

E « laser energy/pulse (joules) 

P . 2 . 

E^ = energy density limit (joules/m ) 

h^ * spaceeraft altitude (ra) 

The second limit to be considered is the maximum receiver aperture. This limit 
results from optical system considerations. The receiver FGV is simply the 
detector diameter (d) divided by the receiver focal length (f). The receiver 
optics F number is the focal length divided by the receiver aperture diameter 
(D„). Thus, for a given optics system F number and detector diameter the 

A 

maximum receiver aperture is. 
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Pjj < d/F (FOV) 


( 2 ) 


The minimum receiver field of view is nominally chosen to be equal to the trans- 
mit beamsidth(of_}. Thus, the maximum receiver aperture can be defined as» 



The radar range equation xdiich determines the return signal energy from the target 
can be written as , 


2 '^-4 ■ 

Ng - Ep n Gif 0 Sj. Sj (4r) R hv (4) 


where, 

Ng ■ signal p-e/pulse 

Ep « transmit energy /pulse (joules) 

Gp * tr^smit antenna gain 

A » irD_/4 (m^) 
e x 2 

a - target cross-section (m ) 

« transmitter optical efficiency 
S. - two-way attoospheric transmissiblllty 
R ■■ slant range (m) 

Ti ■ detector quantum efficiency 
S * receiver optical efficiency 

■ X ■ .._34 ' ■ ■ 

h * Planck's constant (6.6256 x 10 Js) 

■ 14 . 

V * optical frequency (5.66 x 10 Hz @ O.SSpm) 

Thus, substituting 
(5) 


2^4 


For a Gaussian beam, the transmit antexuoia gain is 32/c(^, 
equations (1) and (3) into equation (4) yields, 


K 4 2 j 2 

/ Sv w \ - 


■ MAX * “ ^ ^ « P 
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The design range, R, is chosen as Che range at which the elevation angle of 
Che spacecraft, when viewed from Che target. Is E (nominally 20** minimum). 

I 

/ 2 2 2 

R »-r sin E +*(r +h ) - r cos E < 

e T' e B e 


where. 


r “ earth radius (ra) 


The next step is to determine the required signal energy to meet the ranging 
accuracy goal. Previously, we determined that the variance of the time of 
arrival estimation error, assuming a matched filter or correlator, was 
given by. 




where , 


Ng ■ signal photo-electrons/pulse 
“ background photo-electrons /gate 


A raised cosine pulse shape was assumec, and the pulse width at the half power 
points was T/2. The nominai gate width is T. In order to meet a goal of 10 cm 


Assume a pulscwidth (FWUH) of 


rms range measurement error, ^ "0.667 ns. 

5 ns. Then, 

**s^ / \ 

57 TM 7 ■ \- 7 r- J- U.A 


This equation describes the minimum signal to shot noise ratio (SNR) which 
will satisfy the 10 cm rms accuracy goal with a matched filter detector. 
The effective background energy received is determined by. 


/ lTDj^\ IT 2 

% " ^ V'T^/ ^ A 


where, 

AX " filter bondwidth (A) 

2- o ■ 

N ■ background radiance (w/m -A-G) 


SYNTHESIS AND ANALYSIS GF 

PRECISE SPACEBORNE LASER RANGING SYSTEMS 


REPORT MDC El 729 
FEBURARY 1978 


Then 


, let FOV ■* and substituting aquations (1) and (3)» ve find. 


- T A A 

D 


^ (ff) 


.005 m. 


At 0.532 the sunlit earth background radiance is approximately 
0.017 v/m -A-ster. Assume AA ■ 10 A, T - 10 ns, n “ 0,25, S^. ■ 0.4, d 
F - 2.5. 

N„ - 1118.53 p-e/gate 

D 

Then, solving equation (8) wo find the required Ng is, 

N ■ 118,8 p-e/pulse 

O ' - 

Next, ussxime h^ * 333 (i80 n«i),E " 20°, “0.8 

Then, R " 838 km , and, 

8 2 

N„ “ 4.18 X 10 E. a/E^ 

®M»X * f 


Next define the performance margin as the ratio. 


K • (Ns ) / (No ) " 3-52 X 10 E. a/E^ 

P M\X ®RQD a P 

The maximum safe visible frequency energy density foi^ human exposure la on the 
^7 £ —5 

order of 10 joulea/cm . In order to be conservative, choose E^ " 10 

Joulea/m (23dB safety margin) . 

_4 

M “ 3.52 X 10 o/E 
P P 

Noto that appears in the denominator of this expression, implying that smaller 
energy per pulse lasers have a greater performance margin. This, of course, 
results from choosing the transmit beamwldth as small as possible and the receiver 
aperture as large as pOBSible, In the experiment, we will choose E ns large 
as practical. A value of ■ 50 mj/pulse @ 0.532 pm appears achievable, 
resulting in, . 

M -0.007040 
P 

The targets are aaaumed to be arrays of cube-corners with a maxlinum radar 

7 2 4 

cross-section of 10 m , resulting in 7.04 x 10 (48.5 dB) . A practical 
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systetu, can be configured which is up to SO dU leee efficient than 

the tnaximum perfomancc aystem. Mote f.hat the traixamlt beamwldth and receiver 
diameter of the maximum performance aystem are 0.339 mrad and 5.9 m, respectively. 
A practical receiver diameter for a Shuttle flight is on the order of 15 cm 
(6 inches)* or *^ 32 dB less receiver area than the maximum performance system. 
However, the received background energy is also reduced by 32 dB, which red^iccs 
the required signal energy by 12 dB, for a net reduction of 20 dB compared 
to the maxiittum performance system. 

We conclude that a maximum performance system would collect approximately 30 dB 
more return signal than would be required to meet a 10 cm rms accuracy goal. 
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2.1,2 Multiple Target Ranging 

The purpose of the laser ranging experiment is to obtain data which will allow 
detetminatlott of the relative distances between targets in a grid. The simplest 
approach, conceptually, is to emit a laser pulse which Illuminates all of the 
targets in the grid, and measure the difference in the time of arrival of the 
reflections from each target. Then, after a first order correction for the 
range rate^ the relative distances between the targets in the direction of the 
observation cat\ be estimated with considerable accuracy. This process Is then 
repeated for a number of observation angles to obtain a three dimensitmal rep- 
resentation of the target grid spacing. 


This approach is reasonable where the physical slae of the target grid Is small, 


but quite difficttit to implement where the grid siae is large. Let be the 
diameter of the minimum circle which juat contains the target grid. Then, the 


minimum receiver FOV Is simply D„/h^, 

<vvV^ 


The optimum transmit beawwidth is 


As discussed in the preceding section, the maximum performance system had 
approximately SO dB of margin above the 10 cm rms performahee goal for a 
SO mj/pulse, 0.339 mrad system with a 5,9 m dia. collector. This would result 
in an illuminated spot on the earth of *v 0^08 km diameter for a 333 km orbit 
altitude. If the transmit beamwldth and receiver POV are Increased to view a 
one km diamter grid, the perfotmance loss is 44 dB, or nearly all of the 
available 50 dB performance margin, in order to rjutge from all targets In « 
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10 km diameter grid, another 40 dB of performance loss could be encountered, 
requiring a total of 4 orders of magnitude Increase In target cross-section 
and/or transmitter power. We concluded tha^ this concept does not appear to 
be a feasible alternative. However, we elected to perform a more detailed 
parametric analysis to determine the limits for usable grid sizes. 

One additional parameter, AA, the optical filter bandwidth, was Included as 
a dependent variable. Although It Is plausible to construct optical filters of 
nearly arbitrary bandwidths for narrow FOV applications, when the receiver FOV 
is On the order of a degree or more, existing filter technology limits the 
minimum usable bandwidth. Figure 1 shows one such projection, which is con- 
sidered typical of the current state of Che art. The smooth curve Is the result 
of least square fitting the available data, and results in 4\ 'v 3 + *155 ot,,, 
where AA is the filter bandwidth (A) and is the filter field of view (degrees), 
there Is a current resurgence of interest in narrow bandwidth optical filters 
with wide usable cone angles. Techniques such as mosaics, selective gas absorp- 
tion, and fiber optics channels are plausible. However, these concepts are not 
yet In the laboratory demonstration phase, thus it was Judged inappropriate to 
demand better than currently available optical filter performance. We assumed, 
then, that the filter bandwidth followed the square law expression, above, and 
that the filter was located, in effect, at the receiver aperture. Note that 
the effective filter field of view is therefore equal to the receiver field 
of view; if a smaller diameter filter were used, elsewhere in the optics train, 
the filter FOV would be greater than the receiver FOV. 

Typical high speed photo-multipliers, such as a static crossfleld PMI, have a 
small (4-6 mm) photo-cathode Submerged Inside the tube, with a modest sized 
window ('^'.5" diameter). Typically, these detectors require the optics system 
to have F numbers on the order of 2.5 or more. 

If the tube were redesigned, it is possible to hypothesize somewhat larger 
photo cathodes with a lower allowable F number, but order of magnitude Improve- 
ments are considered virtually unattainable in e high speed device. For fea- 
sibility evaluation purposes, we chose to aseune a detector diameter of 10 mm, 
and a minimum F number of 1.2. these limits are used at 0.53um wavelength. 
However, at 1.06pm, the situation is somewhat different, since the most attrac- 
tive detectors ere photo-diodes, with or without avalanche gain. These devices 
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are limited to aisee ou the order of 1 am\ diiunoter (active area) in order to 
obtain aut'friclent bandwidth, however^ it la poeeible to uae special techniques 
to decrease the allowable F luimber to mO.S. 

Finally, wo elected to tUitonnlno the required signal energy lor an optimal time 
of arrival estimator rather' than the matched t'ilter used in the previous section. 
For a maximum likelihood (TH) CLtmo*'of-arriV'al. estimator (w'ith a raised cosine 
puise) , the variance of the estimation error iSj, 



This can bo Inverted to obtain, 

2ir/ 

The reqnired signal energy for a ML astl»nitor is only slightly lass ttian for a 
correlatov when the background level Is large. The ML advantiige ixicreases as 


Ng “ K where K « 1^ 
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the background level is decreased, and approaches 3 dB at zero background level. 

For feasibility evaluation purposes, we would like to detennine the performance 
parametrically with respect to the target grid dimensions. From a practical 
standpoint, the transmit energy per pulse is limited to about 50 raj, hence, the 
only unconstrained variable remaining is the target cross-section. Thus, we 
elected to solve for the required target cross-section as a function of grid 
diameter. 

The reference system is based on the guidelines in the SOW, l.e., 50 mj per pulse 
§ 0.53ym wavelength, with a 5 ns (FWHM) pulse width. We assumed the receiver 
aperture was no greater than 0,15 m (6 in.), and a 25 % quantum efficiency 
detector would be used. The transmit and receive optics were assumed to have 
80^ and 402 transralsslbility, respectively. Atmospheric transmisslblllty is 
generally expressed as exp{-K/sin E), where E is the elevation angle of the 

propagation path, measured at the target, and K is a constant on the order of 
.08 (one-way) at 0.53ym wavelength on a clear day. 

Figure 2 shows the results for the nominal system, at spacecraft altitudes from 

200 km to 400 km In 50 km steps, for 20" minimum elevation angle, assuming 

Sunlit Earth background. A reasonable maximum target cross-section is probably 
7 2 

on the order of 10 ra , thus we see that it is virtually impossible to obtain 
satisfactory operation at grid sizes greater than 1 km. Note that these results 
are for zero-margin with an ideal detector (ML). If we allow 6 dB margin for 
real detection processes/ the design point would be for o “ 2.5 x 10 m , which 
limits the grid to sizes between 0.4 km and 0.6 km. 

Figure 3 shows the results, for the same conditions, for the l.ObMm variant of 

the reference system. As usually found, the performance at 1.06pm is worse 

than at 0,53pm, thus little reason can be found to pursue this alternative. 

The major problem with the broad beam approach, of course. Is the reduction of 

energy density at the Earth surface as the transmit beamwldth is increased. 

Thus, we hypothesized a multiple tr^smit beam, single receiver system concept, 

which could, effectively, accomplish the multitarget ranging objective. Figure 4 

shows the results under the 20" minimum elevation angle constraint. Although 
7 2 . 

the 10 m (zero marglh) performance has been extended to 14 to 24 km, signifi- 
cant improvement does not appear plausible^ If the minimum required elevation 
angle is increased to 30**, operation with grid sizes up to 'v30 km appears feasible, 
as seen in Figure 5, but larger grids appear unreachable. 
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Figiire 4 


Fig«re 5 
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One possible alternative Is to design the system for operation only at night. 
Although operationally undesirable, the Improvement in perfonaanee capabilities 
may Justify the decision. Figure 6 shows the results for the wide beam concept, 
with the background radiance set at zero. Clearly, the broad beam concept is 
simply not a viable alternative for grid sires greater than 1 km, even with 
zero background. Figure 7 shows the results for the narrow beam wide FOV concept. 
The influence of the detector size on the receive aperture is clearly seen since 
the required target cross-section is constant until increasing POV requirements 
begin to limit the maximum receive aperture size. 

The concept of zero background level, however, la never really encountered. 

Thus, before we conclude that the concept Is feasible, we must make an 
assessment of the impact of some modest background level on the performance. 

Since the Moon is the brightest object in the night time sky, we can estimate 
the background radiance of the Earth at night due to moonlight as the ratio 


TARGET CROSS-SECTION PARAMETERS 



REPORT MDC El 729 

VOLUME II 

FEBURARY 1978 


SYmHEStS AKO AHALYSIS Of 

PRECISE SPACEBORHE LASER RARfilRG SYSTEMS 


REPORT MDC El 729 
EERURARY V97S 


VOUiMF n 



ot suftcteal V}idt«nG<i tUw r«diH\vc.« ol tU« avmlit RAVtU, 

A . i * 

i,e,> (A.7 X 10“ /28‘ioy x *0l7 ^ x 10 v/m 

Xw ovOtfV to bo oov\i4t>vviVt tv%'» wo Irtoroiiaod x Uo brtObjtVvHmd level tbtoe ovdoi'a 
of nuij^uVtwdfcS wixb tbo toaulta aX\ow\^ io FljSvwo 8. Vfe cowc luvlo tbat nett bet 
tiw wide benttt not^ tbe vtAtvow boinvi* wide VyiV covio^^^ la viable I'ov wodcst to 
liU‘ise §tida. 

VUvoo tit» 0 “*of“atvival ilXtAl ontiwatlon tecbuivieea wove ivtvoatlAatod fov the 
laaov van^tinx oxpoeimont , n mtobod ftltev tocXnvldn^^ nv\d ixiio alldlnii wlnd^^ 
intn^ivitot to.chnlvinoa. The ve levant V'erfon'ivnvco analyaos avo dlacnaaod in tbe. 
next pavajsvivplxn * 

Si piuig ^ ^ ' ' 

A slid tftiix window Intoixtntox^ can bn coitstintotod^ oonon\vtvu!ill>S as shovn in 
Hxttvo 9. litis np\vi‘oaob nnos a tbtoidtold dotoctov to dotnvmlno tbe IMA 


SYNTHESIS AND ANALYSIS OF 

PRECISE SPACEBORNE LASER RANGING SYSTEMS 


REPORT MDC El 729 

VOLUME II 

FEBURARY 1978 


TARGET CROSS-SECTION PARAMETERS 



estimate. An automatic threshold controller (ATC) adjusts the detection 
threshold to minimize £alse alarms. The TOA estimate, for this approach, Is 
not ideal since it tends to be biased If the mean signal level is either greater 
than or less than approximately twice the threshold setting* The concept is 
useful,, however, for analyzing false alarm and missed detection statistics, 
particularly where the range gate is large (compared to the pulsewidth) during 
acquisition. 

SLIDING WINDOW INTEGRATOR BLOCK diagram 



Figure 9 

2-15 




SYNTHESIS AND ANALYSIS OF 

PRECISE SPAGEBORNE LASER RANGING SYSTEMS 


REPORT MDC ET729 
FEBURARY 1978 


VOLUME xr 


During the acquisition inode * the laser ranging system will be searching in 
angular space and time to locate the target* If a false alarm occurs, an 
unnecessary target veriflcatlpn sequence will be executed, resulting in a delay 
of correct target acquisition (after the search is resumed). Thus, a single 
false alarm is simply a nuisance. Only if false alarms are frequent will any 
serious problems be encountered. The nominal acquisition period is on the order 
of 2 to 20 seconds; at ten pulses /secondi there are a maximum of 200 opportunities 
for false alarms. A single trial false alarm probability of l/2O0O will there- 
fore limit the probability of more than one false alarm during acquisition to 
about 1%. 

Previously, we found that the probability of false alarm in a single trial for 
a sliding window integrator could be expressed as^ 




(14) 


where T scan time (sec) 



and 


n^^ * mean hackground rate (p 
T * gate width (sec) 

1 » threshold level (p - e) 


when i»n^ » 

L-1 





e/sec) 



Clearly, for any particular background rate, gate width, and scan time there is 
a threshold setting idilch will limit the false alarm probability to an accepta- 
ble Value. Note, however, that for a spacecraft in low altitude orbit, the mean 
background rate observed with a narrow field Of view detector can vary rather 
rapidly ^ We devised an automatic threshold control (ATC) technique to -cope 
with a rapidly Varying background level . 
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Assume that the signal, if present, will be in a scan time window, During 

a time period, immediately preceding the scan time window, the output of the 
sliding window integrator is fed to a peak detector. The threshold is then 
set just above the maximum integrator output sensed in the seconds time period. 

The threshold setting, thus selected is a random variable. The probability 
of selecting a specific value, is, 

p(L) = e - e ^ - (16) 


Then, we see that the probability of a false alarm during the scan time window, 
assuming is constant during the T^^ + T, time period, is simply, 


PrlFA) 


E a 


e )(e 




“®L-l'^2. 
- e ) 


(17) 


L=1 


This equation has been evaluated numerically. We assumed T = lOr, and varied 

5 ^ 

T^ from lOOx to 10 T for a widerange of background levels (n^T = .001 to 100); 

the results are shown in Figure 10. As can be seen, T.^ = lO^T to lO^x 

—3 

results in a usable false alarm probability (<10 ) over the range of background 

levels examined. The behavior of the false alarm probability curves for small 
background levels reflect the transitions of the pump“Up period as larger 
threshold settings become probable. These transitions become blurred as the 
background level Increases. 

The time available for the pump-up period is the time between departure of the 
transmit pulse and the beginning of the range gate, Which reaches a minimum of 
'^2.2 ms for a direct overhead pass. If we allow "vi ms for the pump-up period, 

5 

and assume t 10 ns, T 2 /T = 10 , and we can expect satisfactory results for a 
100 ns range gate. During initial acquisition, however, the range gate is con- 
siderably larger, more like 3 ps. Figure 11 shows the performance of the 

5 ■ 

ATC for T_ = 300x. Note that for T„ = 10 T, the false alarm rate for ivx “ 10 

I ■ _5 ■ 2 , 

has increased from 6 x 10 to 1.8 x 10 , a factor of 30 increase, showing a 

nearly linear dependence of to both T- and T„ where is small and 

V-^‘ ’■ 

Thus , we conclude that a pump-up ATC with a sliding window integrator is a 
satisfactory technique to Control the false alann probability over wide range 

^ ■ 2-17 
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of bttckgritjund levels. Further, a pump-up time period of 1 ms or less will 
nominally result in a satisfactory false alarm probability * and can be 
accommodated readily. 

The next step Is to coiiq)ute detection probabilities for the pump-up ATC system 
when the return pulse is present. This can be estimated very easily by 
stmpiyconsiderlhg the detect ion probability whore the sigual pulse is con- 
tained tdiolly within the si idle s window integrator. If we ignore false alarms 


‘MD 


OO 

£ 


-a T„ -a *T„ 

(e ^ ^ - e E — - e 


(18) 


K-0 


Tltis equation has been evaluated nuwerc tally with the results shown in Figures 
12 and 13. Comparison shows that increasing has a stwall effect on tho 
required signal energy for a given missed detection probability. Figure 14 
shows the missed detection probability for a stationary threshold system for 
comparison purposes. The threshold was selected to keep f Ji TO and T^ » lOx, 
As can be seen, the difference between the curves for the random threshold with 

■ ji 

T^ * 10 T and the stationary threshold is very small (<2 signal p-e). Thus, we 
conclude that the pump-up ATG system concept is a viable candidate ceebnlq\io 
for the laser ranging system. 


Correlator Performance 

Previously, we determined that a correlator or matched filter technlqnc! yielded 
improved tlme-of -arrival estimates compared to a sliding window Integrator with 
threshold detection. One possible implementation techniqtie is to employ a 
tapped delay line with pulse shape weighted summing of the outputs. 

Figure IS shows a typical example, with six taps. .Since a finite number of taps 
are to be used. It is desirable to determine the degiadation in perforniaivce as 
coitqjared to a continuous matched filter. Consider the case where the pulse 
shape is well approximated as a raised cosine. 

X(t) « tt a + cos 2irt/T), tt|£f/2 (19) 

njr is the moan energy in the received pulse. 
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output' 


OUTPUT TO PEAK 
DETECTOR 


When the received pulse is exactly centered in the delay line, the niean signal 
energy in the jth bin is sinq>ly 

- «g fT/N -I- I sin (tr/N) cos (2 tt - 1/2))] (20) 


where N ■■ niunber of taps. 
Assume the weighing functions are » 


Wj - 1 + b cos (2 tt ( J ~ - 1/2)) 

where b is a constant used to vary the filter characteristics. 


The second moment generating function of the output of the jth tap is, 




W,S 

a J - 1) 


The second moment generating function of the output of the H umme r is. 


^ H.S “ 


J. JU WSn.Tj, WsJL 

* E 'fj • E h.i « ^ ^ E * ^ - E".,: 

j-1 j»l 


where is the background rate in p‘-e/sec. 
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In order to reduce this express ton to a mote compact form, we will use the 
following relationships . 


^Zcose ^ + 2 / ^ I^^(Z) COB ke 


sin jx = sin (ti x /2) sin ((N + 1) x /2)/sin x /2 


S cos jx » sin (N x /2) cos ((N + 1) x /2)/sin x /2 (26) 

W.S 

We first expand e ^ as indicated in the first relationship, and then perform 
the trigonometric summations , and find, 

- nj^TtI^(bs)e® ^ 1] + n^Tl(l^(bs) + I^(bs) | sin tt/N) e® - 1] (27) 

1 11 2 

Then, noting that (0) = n and t|j (0) = a , we see, 

3C 3C 3C. X 

“x = V* V U8) 

®x ’ + b /2) + ^ ^ ^ 

TT N 

Thus, we see that with b = 0, the filter is a simple integrator, and the mean and 
variance are equal, as expected for an ideal phot^^.ectrpn coTmter. Where b “ 1, 
and N is large, > > 

- I V I ‘W 


The next step is to estimate the effect of varying S and b on the performance of 
the system. One of the more significant criteria is detection probability and 
false alarm probability vs threshold setting. 


In order to estimate these probabilities, we use Chernov bounds as given below. 

Pr (X < yJ 5 exp [t|> (s) - * (s)]/y ip„'(s), S < 0 (32) 

3W' A 
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Pr {X > y) < exp (s) - Si|/ (S) - '(S)]/y =■ '(s). S > 0 (33) 

X 3C X X 

We evaluated these bounds numerically, with the results shown In Figures 16 and 
17. Although varying b causes both the false alarm and missed detection curves 
to shift, the effect on detection probability is rather small, as can be seen 
from the intersections of these curves (highlighted with a circle) . We could 
conclude that there Is an optimum b (a function of background level) , but the 
effect Is small enough that optimising b Is not a high payoff actlArlty. 

Varying the numbar of taps does not affect the false alarm statistics, and has 
only a small effect on the missed detection probability as can be seen from 
Figure 17. 

We were concerned, however, that the finite number of taps might cause the output 
of the summer to be somewhat distorted. This was evaluated with a simulation, 
with the results shown in Figures 18 through 20. The delay line was simulated 
with AO discrete time slots, and the effect of various groupings tested with a 
repeatable random sequence. In the first figure, each time slot was weighed 
individually, and the resultant curve of the summer output Is quite smooth. The 
next step was to group the time slots Into 8 taps (five time slots each), with 
a resulting summer output curve shown In Figure 19. Some distortion Is evident, 
but of little significance. Figure 20 shows that for 4 taps, a significant 
distortion Is observed, which could cause problems in determining the peak. 

Thus, we concluded that 8 taps is a reasonable choice. Four additional runs 
were made in the 8 tap configuration using different random sequences, with the 
results shown in Figure 21. Although some roughness is apparent, the pulses are 
reasonably smooth and the peaks detectahle. 

Alternatively, the correlator could be implemented slightly differently such 
that a zero^crosslng detector could be used to detect the pulse center. This 
wduld be accomplished by changing Wj to be, 

Wj - cos (2ir ((i + .5)/N + 0.25) (34) 

Figure 22 shows the summer output as determined by simulation for three pulses. 

The negative transient preceding the zero-crossing would trip a threshold detec- 
tor used to arm the zero-crossing detector. Determination of the performance of 
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the threshold detector Is s bit more laborious in this case* since the pulse Is 
not completely within the delay line, and the high order Bessel functions do 
not reduce to tero. 

Since we previously found that 8 or more taps was essentially equivalent to on 
ideal matehed filter, wo will onnlyre the ‘'continuous" case. The filter output 
where only background is present is describad by the second moment generating 
function. 




cos(2TTa/T) 


li det 


Het\ce, 


(35) 


"b - ® 


Vdien only signal Is present, we find, 


T/2 


\{>s 


(s 


,t) - J n^fl + cos(2lT(t - a)/TUu® ( 36 ) 


t-T/2 


Expanding the exponential In a Bessel function series, and performing the trigono- 
metric integrations results in, 

t<i^(S,t) - u^t| ir (s) - 11(1 - t/T +1^ sin 2irt/T) 

-lj^(8) t(i - t/T) sin 2TTt/T + |(1 - cos 2lTt/T) 1 


*20 


ir 


(>1)^^^(S) 
K(K^ - 1) 




(1 “ cos 2irt/T) 


- (1 - cos 2irKt/T)] sin Kir/2 

« 

t (sin 2irKt/T - K sin 2irt/T] cos Ktt/2J 


(37) 
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- (f + i-) tl„(8) - 1] - + i 

"slV£ V 

XStt “ 5tt 
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- T/4 (G V 284832 t is closer). 

n H t^(a) 

)n^TI^(s)^ 

+ (38) 


Thus. 


“a " ? ^ 

S £ H TT S 


“s" - <? * fej V 


(t - T/4) 




r, 2 3 „ T 

0 T" n T 

8 - 2 8 


Then, the signal to noise ratios for the two filter configurations (for threshold 
detection) are, 

SNR (MF) - (3/2 n T)^/ (5/2 n T) - 0.9 n T 
s s s 

3 1 2 o 1 

SNR (ZC) - t(f + |jp) /(f + 0.5931 n^T 

Thus , we conclude that threshold detection with the zero-crossing detector con- 
flguration requires 1.8 dB more signal then for the matched filter configura- 
tion for comparable performance. The zero-crossing detector Is easier to 
implement reliably than a peak detector; the 1. 8 dB penalty can be accepted 
since we esqiect to normally operate in the large signal region where detection 
probabilities are not a problem. 

Split Gate Sliding Window Integrator 

An alternate concept, to evade precision wide bandwidth tapped delay lines, is 
to use a split gate sliding window integrator, as shown in Figure 23. This 
device Is very similar to the tapped delay line with zeri -crossing detector for 
time of arrival estimation. 








f>-^i 

8 ! 

^1^1 
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SPLIT GATE SLIDING WINDOW INTEGRATOR BLOCK DIAGRAM 
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The false alarm and missed detection probabilities are given by the equations 
derived for the sliding window Integrator previously. The performance for time 
of arrival estimation can be shown to be 
— ^ »2 


T 

16(SNR) 

2 


where SNR ■ s 

T =* total gate width (2*FWHM) 

n. “ background rate (p-e/sec) 

D 

N ■ signal p-e/pulse 

S ■ ■ ■ ■ ■ ' 

and a raised cosine pulse shape was assumed. 

For comparison purposes, we note that If a correlator were used, 

^ 

ot 2 

“ 2ir‘‘(SNR) 

If a maximum likelihood estimator were used, 

— -s- -,2 


ATT [N^ + Oj^T - 2Ng%T + 
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At zero background ratOt these become* 


a 


16N 


S«’ 


COR ”9 


ML 


Under background limited conditions (n^T»N^) , 

7T 


T^(njjT) 

(n^jT) 

T‘(n. T) 
D 

16 N ^ 
8 

* 2 2 
sc 

• 2 2 
eoR 


ML 


Thus* the rma TOA error for the split gate tracker is '^'1.11 times the rms error 
for a correlator operating on the same signal. The ML estimator is either 
better than or equal to the correlator, but not by much u:less the background 
rate is very small. Alternatively, we can express the performance in terms of 
required signal or SNR. In this context, the split gate tracker requires 
'^0.9 dB more signal to shot noise ratio than the correlator for equivalent 
performance . 


2.1.4 Ground Return Pulse 

The emitted laser pulse can be reflected from clouds, the Intended target, and 
the ground. In order to determine the relative magnitude of these reflected 
signals, first consider the signal reflected from the ground. In many cases, 
the ground is well represented as a diffuse (Lambert) scatterer. The magnitude 
of the return pulse is then given by, 

Nq - Pt^SrA^ (rcos^J^R hv (^2) 

where 

S_ ■ scattered signal p-e/pulse 
G 

P^ » transmit energy per pulse (joules) 

- transmittance of transmit path 

S “ transmittance of receive path 

• • R ■ , - . .2 

A • receiver effective area (m ) 

e ■ ■ 

r ■reflectance of the ground 
4 ■ zenith angle 

R ■ range (m) 
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h » Plmnk*8 coMtant 
V •> optical frequflincy 

I£ the aurface Is a more nearly specular scatter at near normal incidence (^•Q)« 

»S * “g * 
where \ 

6^ • full planar heamwidth of the scatterer 

For example, even moderately choppy water has a significant gain compared to 
a diffuse reflector; at near nomal incidence* On the other hand, a concrete 
road is a nearly lambertaln reflector. 

The reflectance coefficient, X, varies considerably* A typical lower limit for 
black earth is *'' 0.02. An upper limit is M)*8, for newly fallen snow. 

At normal incidence (looking straight down), the energy in the reflected pulse 
from new snow is ‘vl4d p**e/pulae, which is easily detectable, although considerably 
smaller than the expected signal from the target at this condition. Thus, a 
potential false acqulaitlon problem exists if wa must reacquire near cenith, or 
if the target is snow covered. Hear the horizon, the return la much amaller 
(■vl.d p<-e/pulse) and grossly spread In time, and is, therefore, relatively 
unimportant. 

In order to estimate where the false acquisition possibility ceases to be a 
problem, consider a laser pulse temporal and spatial envalope given by, 

S|{t,r) » ' ( 43 ) 

The reflected pulse is then, 

S_(t) * / / a exp {— 

0 0 ro 
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(H^tan ♦ ’cos e) 




) de rdr 


(44) 


REPORT MDC E1729 

VOLUME II 

FEBURARY 1978 


SYNTHESIS AND ANALYSIS OF 

PRECISE SPACEBORNE LASER RANGING SYSTEMS 


REPORT HDC El 729 
FEBURARY 1978 


VOLUME II 


• we see. 

Su{0) • / a exp - | de rdr 

^ 0 0 CV 


(45) 


IfitercKanging the order of integration yields, 

2it 


S-(Q) » / (|) — — - - y - 

o ■ / 1 . 4tan A COS 9 

— irr— 


T 


) 


an 


/ t ^ 4 

V *^0 ’ 


5^(0) = 


a n r 


0 


1 + 


4r^ tair 4 


FT 


CN 


the total reflected energy is determined to be. 


» / S„(t)df = dTr- n 


3/Z 


(46) 


(47) 


Thus. 


Sr(0) = N6(4)/T/n(l + (ZTq tan V 


(48) 


Now. S^(0) is the peak power in the reflected pulse. At normal incidence, all 
of the reflected pulse would be detectable in a gate width. However, as the 
nenith angle increases, the pulse is stretched in time, and the maximum energy 
in the gate is approximately S^(0)T. Thus, we can compute N^(4) from Equation (48), 
and then compute S^(0)T to estimate the peak detected signal in the sliding 
window integrator. 

This will yield a good (but conservative) approxinuitiQn if T. the gate width, 
is ndt too much greater than 2x, »ad the zenith angle is large enough to have 
caused a significant decrease in peak intensity. For a 4ns FWIM pulse. 
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T ■ 2,4 ns. Then, for T • 8 ns. 

hs ■ 333 km, 

FOy “ .5 mrad, r » 0. 

estimates of S^(0)T shown in Table 1. 



TABLE 1 


SIOIAL BACKSCAXTER VS ZENITH ANGLE 

ZENITH ANGLE 

N(j(*) 

Sr(0)T 

(DEGREES) 

(P-E/PULSE) 

(P«E/GAIE) 

0 

148.0 

148.0 

1 

147.9 

66.9 

2 

147.7 

34.1 

5 

145.9 

13.5 

10 

139.7 

6.35 

20 

117.05 

2.47 

30 

85.67 

1.05 
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Thus, we see that the energy In a gate width falls off quite rapidly, and becomes 
virtually Insignificant as the zenith angle Increases above 10*. At near normal 
Incidence, we could also experience a rise In detected background due to specular 


scattering. However, the likelihood of conducting ranging measurements at very 
low zenith angles Is quite small at best and can be completely evaded by 
operational constraints (l.e., look at other targets). 

It does appear possible to use the* system In an altimeter mode If such Is desired, 
although It would probably be necessary to alter the time of arrival decision 
strategy^ and the transmit be amwldth and receiver field of view for this mode. 

Also, since most terrains are about an order of magnitude less reflective than 
fresh snow, It would probably be hecessary to increase the receiver telescope 
diameter to obtain fairly confident ground return detection. Finally, we note 
that the mean spacing between Illuminated spots for a 10 pps laeer Is about 
740m, compared to a beam diameter of <^'165m, thus successive pulses are terrain 
Independent and the range gate would have to be wide enough to cover surface 
roughness. 

In order to test this hypothesis, we estimate the advantage of signal power to 
solar background power* This is readily accomplished by noting that the solar 

■■ O' ' * '■■■ 

spectral irradlance In the vicinity of .53 urn is about 0.2 w/m -A. If we use 
a 5A optical filter, and assume 0.5 mrad field of view from a 333 km (180 nml) 
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spacecraft^ the incident solar power on the ground that is within our field 
of view and optical bandwidth is ^21800 watts* For a 4 ns FHHM» 10 inj laser 
pulse* the peak power Is 2.5 x 10^ watts* thus the peak signal power to mean 
background power ratio is '^115 to 1* and the signal should be detectable above 


the background in a modestly sized gate width. 

Thus* we conclude that ground reflected signal need not be a serious concern if 
we limit the observation region to elevation angles less than 80 to 85*^ for 
normal terrain in the vicinity of the target. Further, it may be possible to 
configure the experiment to allow operation as an altimeter* when not raging* 
at least over favorable terrain. 


2.1.S Laser Ranging Experiment Parameter Selection 

Two types of lasers were considered for the laser ranging experinent* referred 
to as long pulse and short pulse lasers. The long pulse laser is a Q-swltched, 
ca>d^ty dumped NdiYAG laser with an output pulsewidth on the order of 4 to 6 ns 
(FWHM). The short pulse laser is a Q'-switched, mode-locked, cavity dumped 
laser with an output pulsewidth on the order of 0.1 ns (FWIIM). The experiment 
design approach for these two lasers will differ considerably. The longer 
pulse laser system design will be driven by the need to obtain 2 to 10 cm rms 
range measurement accuracy. The short pulse laser system design is almost 
totally driven by detection statistics; if the return pulse is detected* the 
theoretical error will be less than 2 cm rms. The two concepts are therefore 


separately discussed. 

Long Pulse Laser Ranging Experiment 

In the preceding sactions the performance of the laser ranging experiment* 
expressed in terms of theoretical ranging measurement accuracy* was found to 
be a function of the signal to shot noise ratio. 


SNR - 


Nc 





Where* 


2 «4 


Ng - o (E G^) (n Sj S^/(4ir) IT hv) 


(50) 
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Equations (4) and (9) from Section 2.1.1 have been rearranged In (49) and (50) 
to group transmitters and receiver variables, and the more or less fixed 
parameters* 

The transmitter variables are the energy per pulse (E^) and the transmit antenna 
gaiii (6^), vhlch is a function of the transmit beam width (32/a^^). The receiver 
variables are the receiver aperture effective area (A ), the optical filter 
bandwidth (AX), and the receiver field of view (fl^ ■ |- (FOV)^). The target 
radar cross-section (a) is also a variable for this study. The remaining 
parameters are largely determined by the state-of-the-art and the environment. 

The transmit beamwidth and the receiver field of view are primarily constrained 

by the achievable pointing accuracy. Let p be the maximum design transmitter 

-2 

pointing error and let be the transmit beamwidth at the e pwer^ points . 

il 

This would be maximized if * pi 


The wealcest design signal detected would be proportional to e . 




The transmitter pointing error consists of a random component, which can be 
considered as a pulse to pulse variable, and a "static" component, which IS 
relatively constant during a pass. If we assume the random component consists 
of two orthogonal, normal, random variables With equal variances, there the 
probability of the pointing error exceeding p is given by. 


{>p} - Q 




where n " Static podjiting error 


and Q (a, 3) > Harhum Q function* 

<o»e) a 


09 


-( a ^+$^)/2 


K - 0 VB 


The Harkum Q function is 


easy to evaluate numerically, but we seek 


an insight into the relationship between n, p, and o^. For a fixed value of 
Q Cn,3) V 0.1 (one of every 10 pulses is weaker than predicted), a reasonable 
approximation Is , 
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where 


p “ ^ 

K “ 1-12 

O 


(52) 


For Q (o,0) ■ O.OI» the epproxlmatlon is the eame foxia but « 2.0. For 
Q CctfS) ^ 0.001» use "•2.72. Assume - 0.05 mrad, n •" 0.09 sir ad. Thea, 
a value of p * 0.159 nrad will result In ^ of 10 pulses at nooiiiial or greater 
signal strength, and 0.45 ntrad would be optimal. For 99 of 100 pulses 
at nominal, p ■ 0.203 Sirad, and a,^ • 0.574 mrad would be optimal' A value of 
>■ 0* 5 mrad Is a reasonable con^rpmlse. 

The receiver field of view needs to be somewhat greater than twice the selected 
value of p, since there is some net pointing error between the receiver and 
the transmitter. For preliminary design purposes it is, therefore, reasonable 
to choose the receiver field of view equal tO the transmit beamwidth. 

The ''fixed" link parameters Include quantum efficiency, the transmit and receive 
optical efficiencies i atmospheric transmlssability, the range, and the back- 
ground radiance. At 0.532um, a 25% quantum efficiency is achievable in a high 
speed photomultiplier. The values of “ 0.8 and ■ 0.4 have been assigned 
to allow some margin for degradation. The atnuspheric transmissability is a 
function of the optical thickness fT^) of the atmosphere and the elevation angle. 
For targets at or near sea level, an optical thickness of 0.425 at 0.532ym 
corresponds to clear day conditions. Then, S * eXp I-iT^/slnE}* where E is 
the elevation angle. The Shuttle altitude is assumed to be 333 km (180 nmi), 

■2 “ ■ 

and the backgroiiflad radiance, 0.017 w/m ^A-ste correeponds to sna^llt earth 
viewed from space. The design minimum elevation angle was chosen to be 20*’, 
resulting in p 0,0833 and R « 838 km. 

The remaining variables, E^, F, A^, and AX can be choseh to yield any desired 
theoretical ranging accuracy requirement. Assume a 10 cm rms accuracy require'' 
ment, using the split gate slidtog window integrator technique. The requited 
SNR <•■14 (11.4 6dB). Figure 24 shows the required receiver aperture diameter 
to achieve this SNR as a function of the product of energy per pulse and target 
radar eroaS''8ectlon for various optical filter bandwidths. Note thAt t^ 
receiver diameter is inversely proportional to the required rms ranging accu- 
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racy • thus to achieve 2 cn rins the receiver diameter would have to be 5 times 
greater than shown in Figure 2A. An upper limit of 19 cm diameter (7*5 inches) 
was cohaldered reasonable for the receiver* As mentioned previously* a SO mj /pulse 

. 7 2 

laser is considered achievable* and a target radar ccoss-section of 10 m is 
considered feasible* These limits are also shown in Figure 24* to emphasize 
the available design margin. 


REQUmEO RECEIVER DIAMETER FOR 10 CM RANGING 
AT 200 elevation ANGLE ’ LONG PULSE LASER 



103 104 105 106 

Ep o (JOULES - m2) 


in order to meet the minimum acceptable performance* 10 cm rms with a 19 cm 

4 2 

diameter receiver aperture* the product must be at least 1.5 x 10 joules-m , 

?;■■■■■'■ y 2 

asauming a one Angstrom optical filter bandwidth. If e » 10 m * a 1.5 mj/ pulse 
laser would Suffice (zero margin). Alternatively* for a 50 mj /pulse laser 


and o V 10 m ^ a 5.5 cm (1.4 Inch) diameter receiver aperture would yield 
10 cm rms (zero margin) ranging accurac^^^ Note that the required E cr is affected 
by the mlnimtm elevation an8)« For exatople* increasing the minimum 

elevation angle to 30* would reduce the required product by a factor of 
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7*2 (3>6 dB). Decreasing the mlnimuta elevation angle to 10“ would increase 
the required E^a product by a factor of 56.1 (17.5 dB). 

The Shuttle altitude also affects the required E^o product. At 20“ minimum 
elevation angle, increasing the Shuttle altitude to 500 Ku would increase the 
required product by a factor of 4.1 (6.1 dB). 

The selected design values, 19 cm diameter receiver aperture and 5 x 10^ 
joiiles^^, provide greater than 10 dB margin over the minimum requirements, 
which is considered a comfortable margin for preliminary design. 

Table 2 summarizes the link margin for the selected system parameters. 

Short pulse laser ranging experiment - The major concern for the short pulse 
laser ranging experiment is the detectipn statistics; if the pulse is detected, 
the theoretical time of arrival estimation error is well within the most 
ambitious goal. We selected a 0.5 x 10 false ^arm probability goal, and 
a 99Z detection probability goal, and parametrically evaluated the effects of 
optical filter bandwidth and the energy/pulse*^target cross-section product 
on the required receiver aperture diameter, as shown in Figure 25. The 
results are quite similar to the results for the long pulse laser ranging 
racperiment discussed previously, except that the optical filter bandwidth 
effects are apparent even with relatively large energy-cross section products. 
Table summarizes the link margin analysis for a specific point design 
similar to the design for the long pulse laser ranging experiment. 

2.2 Mission Simulation 

The question pf acquisition and tracking cannot be treated in isolation. 
E^eriment characteristics, especially the statistics of the errors in results, 
is affected by selection of the pointing system. In the design phase, the 
experiment behavior can be used to assist selection of fundamental pointing 
system parameters. Among these are Sample rates and allowable acquisitipn 
times. A computet simulation of the experiment is necessary to the design of 
the laser pointing system. 

The experiment has an estimation procedure associated with It. This in fact, 
is the very essence of the ekperlment. Because of the complexity of the System, 
the large number of corner raf lectors which may be considered and the extreme 
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TABLE 2 


LINK MARGIN ANALYSIS 
LONG PULSE LASER 


PARAMETER 

VALUE 

COMMENTS 

1. Transmit Energy 

-13.01 dBJ 

50 mj /pulse 

2. Transmit Lasses 

-0.97 dB 

802 Transmittance 

3k Transmit Antenna Gain 

81.07 dB. 

0.5 mrad beanwidth 

4. Pointing Loss 

-4.34 dB 


3. Free Space Loss 

-265.93 dB 

838 KM 

6. Atmospheric Loss (Dosn) 

-5.40 dB 

0.425 Optical Thickness 

7. Target Gain Product 

206.47 dB 

rn7 2 
10 m 

8. Atmospheric Loss (Up) 

-5.40 dB 


9. Free Space Loss 

-265,93 dB 


10k Receive Antenna Gain 

121.00 dB 

19 cm Dla. 

11. Receiver Losses 

-3.98 dB 

40% Transmittance 

12. Received Signal Energy 

-156.42 dBJ 


13. Energy/Photon 

-184.26 dBJ 

0.532 |im 

14. Received Photons/Pulse 

27.84 dB 


15. Quantum Efficiency 

-6.02 dB 

25% 

16 . Received P-E/Pulse 

21.82 dB 

152. p-e/Pulae 

17. Background Radiance 

-17.70 dBW 

0.017 W/M^-A Ster. 

18; Receiver FOV (Sterrad. ) 

-67.07 dB 

0.5 mrad 

19. Optical Filter Bandwidth 

6.99 dB 

5 A 

20. Receive Antenna Area 

-15.47 dBM^ 


21. Receiver Losses 

-3.98 dh 


22. Received Background Power 

-97.23 dBW 


23. Energy /Photon 

-184.26 dBJ 


24. Received Photons/Sec. 

87.03 dB 


25. Quantum Efficiency 

-6.02 dB 

1.26 X 10 p-e/sec 

26* Received P-^E/Sec^^^ 

81.01 dB 

27. Receiver Gate Width 

-80.00 dB 

10 Nanoseconds 

28. Background P-E/6ate 

1.01 dB 

1.26 p-e/Gate 

29. Received P-E/ Pulse 

21.82 dB 

152 p-e/Pulse 

30. Rec(uired P-E/Pulse 

11.83 dB 

10 cm nos 

31. Link Margin 

9.99 dB 
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TABLE 3 

LINK MARGIN ANALYSIS 
SHORT PULSE LASER 


parameter 

VALUE 

COMMENTS 

1 . Tranamlt Energy 

-16.99 dBJ 

20 mj /pulse 

2. Trensinit Losses 

-0.97 dB 

80% Transmittance 

3. Transmit Antenna Gain 

81.07 dB 

0.5 mrad beamwidth 

4 . Pointing Loss 

-4.34 dB 


5. Free Space Loss 

-265.93 dB 

838 KM 

6. Atmospheric Loss (Down) 

-5.40 dB 

0.425 Optical Thickness 

7. Target Gain Product 

206.47 dB 

1 ft7 2 

10 m 

8. Atmospheric Loss (Up) 

-5.40 dB 


9 . Free Space Loss 

-265.93 dB 


10. Receive Antenna Gain 

121.00 dB 

19 cm Dla. 

11. Receiver Losses 

-3.98 dB 

40% Transmittance 

12. Received Signal Energy 

-160.40 dBJ 


13» Ehergy/Photon 

-184^26 dBJ 

0.532 tan 

14 . Received Photons/Pulse 

23.86 dB 


13 • Quantum Efficiency 

-6.02 dB 

25% 

16» Received P-E/Pulse 

17.84 dB 

60.8 p-e/pulse 

17. Background Radiance 

-17. 70 dBW 

0.017 W/M^ A-Ster. 

18. Receiver FOV (Sterrad) 

-67.07 dB 

0.5 mrad 

19. Optical Filter Bandwidth 

6.99 dB 

5 A 

20. Receive Antenna Area 

-15.47 dBM^ 


21. Receiver Losses 

-3.98 dfi 


22. Received Background Power 

-97.23 dBW 


23. Energy /Photon 

-184.26 dBJ 


24. Received Photons/Sec. 

87.03 dB 


25. Quantum Efficiency 

-6,02 dB 


26. Received P-£/Sec 

81.01 dB 


27. Receiver Gate Width 

-96.99 dB 

0. 2 ns 

28. Background P-E/Gate 

-15.98 dB 

0.025 p-e/Gate 

29. Received P-E/Pulse 

17.84 dB 

60.8 p-e /Pulse 

30. Required P-E/Pulse 

10,93 dB 

12.4 p-e /Pulse 

31. Link Margin 

6.91 dB 
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REQUIRED RECEIVER DIAMETER FOR 99% DETECTION PROBABILITY 

IMO 
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s 

I 

c 

o 
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1 

103 10^ , 10® to® 

Ep (T (JOULES ^ m2) 

Figur* 25 



accuracies desired > the estimation procedure is a very Involved and complicated 
program. Mathematical simulations of this exist but are prohibitively expen^- 
slve in computer time and In cost to be used for engineering design. Because 
of this a simplified simulation was conceived so chat almpst any effect can 
be Included in the statistics If necessary even though the dynamics of the 
nominal trajectory cannot be used for prediction purposes. This simulation has 
been Implemented on a 6500 Cyber computer. 

2,2.1 Program Description 

The basic kinematics of the simulation represents a spacecraft travelling in a 
circular orbit with specified altitude, node. Inclination and anomaly from the 
node at time aerp. This was implemehted by assuming perfectly circular motion 
with uniform velocity. This was to avoid any necessity of numerical Integration 
In the simulation^ The array of reflectors on the Earth's surface were simulated 
by Circular motion representing the earth's rotation. The position of any 
reflector at any specified time Is obtained by calling a subroutine. In this 
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subroucine the positions of every reflector at time zero are calculated just 
once when the subroutine is first called. The program is so dimensioned that 
any number of reflectors can be considered. Rearrangement of reflector grid 
geometries is effected by simple modification of the subroutine. 

Simulation of the sampling is also handled by means of a subroutine. This 
selects the reflector which is to be observed in accordance with any observa- 
tion policy which one may wish to investigate. The time of the observation is 
also provided in accordance with any sample rate specifications. Acquisition 
delays are taken into consideration by having the routine suppress observations 
for a specified time after the acquisition reflector has arrived within view. 

The reflector grid selected for this study consisted of five reflectors 
situated at the corners of a 25 km square with one reflector at the center 
which is at latitude 33 deg north with the sides of the square situated east- 
west, north-south. The purpose of the experiment is to determine the relative 
position of these reflectors each with respect to another. Four "outrigger" 
reflectors were added. It is not required to determine the positions of these 
with high precision. These were added to investigate their ability to aid in 
spacecraft position determination. These four reflectors were positioned at 
the comers of a 200 km square with either its sides north-south, east-west or 
Its diagonals north-south, east-west. 

With these reflector arrangements the measurement policy adopted was to range 
upon the most remote outrigger which is in view (in view meaning that the line- 
of-slght is greater than 20 degrees in elevation) if any of the target grid 
reflectors are hot in view. If the target reflectors are in view, then the 
range to each of them is measured in succession in batches of five measurements. 
After this the outriggers are observed once each and the process is then 
repeated. The intervals between observations are one second unless the observa- 
tion changes to a hew reflector, then it is five seconds. An acquisition time of 
ten seconds was assumed; Subsequeiit evaluation indicated that these assumptions 
are adequate for the laser ranging system design. 

A gimbal arrangement was assumed. This assumed that the system was aboard a 
space shuttle which had its longitudinal axis directed normal to the orbit plane 
and rolled to an attitude so that at the time of closest approach to the central 
ground reflector i the roll gimbal was at Its neutral position. Computations of 
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glmbal angle time histories was provided In a subroutine. Glmbal rates at 
sample times were also computed as was range, range rate and the magnitudes 
of the corrections In return pulse arrival time due to atmospheric refraction. 
If E is the elevation angle and A « 2.4m, B ■ 0.0025m then the range correction 
is given by 

AR <• A/slh E - B/sln^E 
and the atmospheric time of arrival by 
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and the first order velocity correction by 

2 £ A 

c c 

where p Is range and c Is the velocity of light. 


2,2.2 Statistics 

In order to describe the statistics of the estimation procedure which consti** 
tutes the experiment we must specify the state vector which includes all of the 
variables pertinent to the problem. We include the three position components 
of each reflector and the position and velocity of the spacecraft. If there 
are M reflectors then we have specified 3M + 6 variables so far. To these we 
shall add five geopotential terms. According to our assumption of a purely 
circular orbit these geopotentials are considered to be nominally zero. Their 
purpose is only to degrade the statistics of the spacecraft state. The range 
measurements are considered to have an error which Is nearly systematic. An 
autocorrelation function of the form 


R(t) 





(53) 


was assumed with T 20 seconds. This can represent attitude fluctuation of 
the spacecraft. The range is also assumed to have a random error which is 
uncorrelated between samples. The nearly systematic error is appended to the 
state vector. Also appended is a purely systematic error in the longitudinal 
component of the centered mass of the spacecraft. Thus, the total number of 
state components is now 3M -(- IS. 
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2,2,3 Initial Covariance 

The Initial covariance for the target location will depend upon the surveying 

accuracies with respect to the dynamical center of the spacecraft orbit. We 

assume that the center of the grid will be known only to within a spherical 

distribution with a standard deviation of a “ 100 ft. Relative to this bench 

o 

mark the Individual reflectors were assumed to have a circular distribution 
which is uncorrelated between reflectors with standard deviation “ ,3 ft. 
Thus, the Initial covariance for the reflectors is of the form 
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The Initial covariance of the spacecraft state which was used corresponds to 
good tracking (l.e., within two orbit periods) by a ground network. This is a 
complete six by six matrix in orbit plane coordinates. This matrix is as 
follows (in ft and ft/sec).* 


*Refererce (1) Lyndon B, Johnson Space Center Memorandum to FM3/Mlsston Analysis 
Branch from FM8/Mathematical Physics Branch, dated 29 May 1975* 
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.275671+05 -.802681+05 .294705+04 .932597+02 

.799405+06 -.19l430f05 -.79560frH)3 
.116397+05 .182926+02 

.803689+00 

Sjnnmetrtc 
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-.177456+02 

.575508+02 

.128867+01 

-.680565-01 

.240196-01 


-.126507+01 

-.106914+03 

-.334140+01 

.977056-01 

.136106-02 

.525264-01 


The covariance of the spacecraft state Is reinitialized to this at the beginhing 
of every pass over the target grid area. 

The geopotentidl terms considered were G(0,0), S<7»4)» C(6,5)» C(7,6) and S(5»5). 
These were selected by computing the disturbance aoceleratlon acting upon the 
spacecraft by a 10 value Of the expected error in the known values of each of 


the harmonics up to C(7^7), S (7, 7) acting individually. Those selected gave 
the largest acceleration while the spacecraft was closest to the reflector 
grid. The 10 value of the errors in these where considered the apriori 
a tiUidard deviations. These are 


cr “ 0.1 X 10"^ for G(0,0). 

0.5 X 10“^® for S(7,4), 

0.1 X 10^® for C(6,5). 

0.1 X 10“^® for G(7»6), 

0.3 X 10^® for S(5,5) 

Great accuracy in these quantities is not critical to the simulation since these 
are used merely to couple with the spacecraft state components. This is needed 
to prevent a too optimistic spacecraft state coiq>onent determination. 


Initial covariance for the range bias was assumed represented by 
and for the spacecraft center of mass by 0^, ■ .5 ft. 




2.2.4 covariance Propagation 

The initial Covariances which have been discussed must , in order to be used by 
the estiiiuttion procedure* be propagated to the time at which a range measurement 
is taken and also propagated between measurements. For the reflector positions 
this Involves Only the kinematics of the Earth rotation. For the spacecraft 
state this must involve circular orbit state tfahsltlon matrices and the 
coupling from the geopotential effects. The geopotential effects are considered 
stationary . For the range bias the range error autocorrelation f unc tion must 
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be taken Into account. The spacecraft center of tnass error Is stationary. It 
is desired that numerical Integration of Rlccatl equations Is to be avoided. 

Earth rotation Influence upon the reflector covariance Is given simply by 

r 

/ , 'j 

I cos WgAt -sin WgAt “\ i cos 4 )gAt -sin WgAt 0 \ 


sin WgAt 


cos U)„At 


0 


cos o)_At -sin 0 


sin w„At 

It 


cos WgAt 


Where u)_ is the Earth rotation rate. At la the time Interval over which the 
E 

reflector covariance must be propagated. 

The spacecraft state covariance in the absence .of geopotential effects is 
propagated by means of the circular orbit state transition matrix given by 


2-c 

8 

0 

8 ^ 

w 

sO--=> 

0 

-3t»)At + 2s 

2c-l 

0 

-a-c) 

« 

^(-3wAt+4s) 

0 

0 

0 

c 

0 

0 

£ 

u 

-u>(“3U)At+s) 

u)(l-c) 

0 

2-c 

-(-3(i)At+28) 

0 

-H»)<l-C) 

-bJs 

0 

-a 

2c-l 

0 

0 

d 

-Uls 

0 

0 

c 


Hare b) is orbital angular velocity* c « cos wAt, s » sin wAt, and the coordinates 
are respectively attitude* downrange distance* normal to orbit plane distance and 
their time derivatives. The propagated spacecraft state covariance is 


4 , S d> 
o s^o 

where is the initial spacecraft state covariance . 

To include the coupling of the geopotentlal uncertainties a matrix 4*^^ was con- 
sidered. This is given by 

6x5 / 3x5 \ 


( 55 ) 
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here a la the spacecraft acceleration and h represents the geopotential tenns. 
This Is derived in Appendix A. in this case the propagation of both the space- 
craft and geopotential terms is given by 


T 



Covariance propagation for the range error involves the autocorrelation func- 
tion 

R(T) » 

in this case a propagation over time At requires that the covariance of the 

-2 At/T 

range error be modified by first multiplication by e and then increased 

2 -2 At/T 

by addition of the term cf_ (1 - e ). Furthermore every correlation term 

between the range error and the state must be attenuated by multiplication by 
-At/T 

e . One can see that over a long time interval the range error covariance 
is nearly reinitialized* 

The covariance of the spacecraft center of mass error is as mentioned previously » 
not modified during state propagation. 

linmodelled drag effects were also included. These were treated in much the 
same manner as the geopotentlal terms with the exception that no state variables 
were added to the system. Ah unloiown but nominally zero coefficient of an 
acceleration acting opposite to the spacecraft velocity was considered. An addi- 
tional spacecraft state covariance was computed baaed upon the covariance of 
this drag coefficient. This term is given by 

AS - - 4>(t2)0^(t^)ii^(tj^)]Oj^^Iij»(t2) - <t>(t2)r^(tj)4»(tj)l^ (56) 

<idiere i^<t) is now the coupling term expressing the influence of the drag term 
acting from the initialization time of the overfly pass to the time t. The 
tern AS Is the additional covariance due CO unmodelled drag tdiich has been 
accrued In the time interval beginning at and ending at These computa- 
tions have been relegated to a subroutine. 

The estimation procedure which incorporates the measurements and modify the 
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state error covariance was Implemented as a Kalman filter. The calculations 

pertinent to the covariance simulation involve the state covariance S prior 

to the measurement update, the variance of the random error on the (range) 

2 

measurement and the "observability" matrix M. 


M — statCx T 
" ~ 'a range-' 


(57) 


With this nomenclature the updated covariance S is given by the Kalman formula 


s'*’ - s' - s“m^(MS“M^ + MS" (58) 

Some manipulations In matrix algebra shows that this formula Is mathematically 
equivalent to the least squares formula 

(s'*')"^ - (s“)"^ + 

M 


We adopt the former since no matrix inversions are Involved and our application 
does not Involve precise calculation which would be involved In an actual 
estimator. 


The simulation also allows for any component of the state to represent a 
"consider" variable* that Is one which the estimator does not attempt to Improve, 
this can be done simply by changing the covariance of the consider variable 
back to the values which existed before the update. 

2.2.5 Output of Program 

With all the previous calculations* the simulation provides a time history of 
the trajectory with all of the attendant quantities such as glmbal angles 
and rates together with a time history of the entire state covariance matrix. 

Of more direct Interest is the statistics of the error in position of one 
reflector (reflector 1) with respect to another ref lector (reflector j ) . If 
is the 3 X 3 submatrix giving the covariance of ppsltlon components of 
reflector 1 and is the 3 x 3 submatrtx of correlative terms relating 
reflectors 1 and j then the required covariance of position difference is given 
,by.;. . ■ 

®ll * ■ ®1J ■ ®11 
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An average strain error covariance for M reflectors can be fabricated by summing 
for pairs of ref lectors . This is expressed as 

M 

®avg “ ~M(M-l) ^ " ®ij^ 

' ■ ^ ■ U1 ' 

Also obtainable from the state covariance matrix is an estimate of obtainable 
pointing error. The derivations of the computations for this are presented in 
Appendix B. We must note that this pointing error is based upon the covariance 
of the simulation under discussion and not necessarily the covariance expected 
to be provided by an actual on-board mechanization of the navigation accuracies 
derived from ranging. 

2.2.6 Results 

The covariance simulation vhich has been described was applied to a shuttle 
borne ranging device travelling in a 180 mm altitude orbit having a 55 degree 
Inclination. The node was selected so that the first pass-travelled across a 
point 1 degree east of the central reflector of the array which was at 33 
degrees of latitude. The outrigger reflectors were 200 km apart on a square 
with north-south « east-^est sides. The spacecraft was travelling northeastward 
at this point. The node was reinitialized for a second pass over this same 
point but travelling southeastward. The ground tracks for these are shorn in 
Figure 26. the measurement policy previously described was used. Flg\ire 27 
presents the ginbar angle history for the first pass^ This suggests that a 
renumbering of some of the reflectors might reduce the maximum gimbal rate 
slightly but the rate requirements are not serious . Range history is presented 
in Figure 28. 

Presentation of covariance results in terms of 3x3 matrices or just their 
diagonal terms can sometimes obscure the tfi^ nature Of tlie accuracies obtainable 
because Of the inappropriateness of the coordinate system selected. In order to 
better see the degree of accuracy preaent in the reflector relative position 
covariancej the eigenvalues of the coverianCe matrix were generated. Although 
time histories of these as the ranging process continues represents a principle 
axis coordinate system which Is changing with timei the results of doing this 
still provides a better Indication of system performance. 
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RANGE TIME HISTORY 



A curaory examination of the experimental procedure would lead one to the con- 
clusion that the position difference accuracy between a pair of reflectors 
would be on the order of the random error of the ranging measurements and 
that both the range bias and the spacecraft position error should cancel in the 
differencing process. Indeed a one-dlmensional model of the experiment with a 
stationary observer and one range measurement on each reflector will give a 
range difference variance of 


2a„ 




cf 

2 m 


(61) 


where 0 ^ is the a priori variance of reflector position. However, the simula- 
tion did not behave this way with target reflectors having an angular separation 
as viewed from the apacecraft. The range bias cancels but the spacecraft pOsi- 


t ion error does not. Appendix C shows why. In this case the system behaves as 


2 2 2 2 
Vpi. + 2o_ ) sin fi instead of lust cr 

2^ Sy m 

where is the separation angle and o_ if the variance of spacecraft position 

°y 


if the random error variance is 0 4 

m 
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perpendicular Co Che line of sight and In che plane of Che reflectors. Since 

o_ can be quite large » a- sin 5 can easily dominate over a . 

&y by . m 

Furthermore additional range measurements dp not improve this since spacecraft 

position error is systematic during the observation period. 

Figure 29 shows the principal io position difference accuracies for the relative 
positions of reflectors 1 and 2. These are the square roots of the eigenvalues. 
Notice that the minimum eigenvalue does not immediately decrease as the reflec- 
tors in question are viewed but must await some reduction in spacecraft posi- 
tion. The intermediate eigenvalue begins tO' be reduced when sufficient angular 
travel makes another direction observable. The position difference notmal to 
the orbit plane requires a second pass which observes this component. 

The Spacecraft position determination which is involved in the first pass is 
shown in Figure 30. It seems apparent that a spacecraft accuracy of less than 
1 ft can only briefly be obtained near closest approach. This is the limiting 
effect in experiment accuracy. Modification of the sampling policy cannot 
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appreciably Improve this with the possible exception of additional Independent 
tracking of the spacecraft. Any additional orbit deteniiinatlon must measure 
to an accuracy of greater than 1 ft referenced to the reflector grid. This is 
the rationale for the extra "outrigger" reflectors positioned so as to provide 
good triangulation for spacecraft position. 


SPACECRAFT POSITION DETERMINATION 



The effect ef moving the outrigger reflectors closer to the target grid was 
investigated. Placement of these four reflectors at 100,000 meters rather than 
at 200,000 m degrades the strain accuracies in the target grid by about 30%. 

The outrigger positions have been subsequently changed to be at 200,000 m but 
situated north-south and east-west of the target grid. This is to provide more 
lateral viewing of the spacecraft because of its 55* inclined orbit over the 
area. ' . 


The pointing uncertainty due to the uncertainties remaining in the relative 
spacecraft reflector positions is shown in Figure 31. 


■V Ir S . j' .ir 
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ERROR IN KNOWN DIRECTION FROM SPACECRAFT TO TARGET 



The simulation was exercised to determine the changes which would result from 
more remote encounters of the spacecraft with the reflector grid. This was 
done by varying the orbit node so that orbit ground tracks would have specified 
longitude differences vri.th the reflector grid of the grid latltutde. This was 
done mainly to check ranges and glmbal angle histories to be expected. Figure 32 
shows that the accuracy deteriorates rapidly after 4 degrees of separation. 

Figure 33 shows the giribal rate histories for a longitude separation of 6 degrees. 

Further mission analysis studies are not within the jurisdiction of the contract. 
All that was needed was corroboration of the selection pointing system features 
from ah experiment accuracy point of view. The results so far indicate that 
the pointing system design is compatible with experiment objectives. 

2,3 Acquisition Process Analysis 

The expected angular uncertainties due to ephemerls and target location errors 
significantly exceed the pointing uncertainties due to attitude reference and 
pointing control errors, at least until the first Kalman filter update is accom- 
plished* The acquisition process consists of a sequential search over the region 
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Figurt 32 
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of uncertainty until a valid target la detected* At that point, the angnlar un- 
certainties are reduced, and the normal ranging mode commences* 

Let fo(X) be the aprlorl probability distribution of target location, where X 
is a two dimensional angular vector. Let g (X - X^) be the probability of de- 
tection, given the target is at X and the beam center is at X^ (for nth trial). 
Then the probability of detection ^ the nth trial, , is, 

P (n) - 1 - / fo(X) n U-g(X-X^)I dx (62) 

1"1 

The probability of detection ^ the nth trial, Pd (^) * simply , 

n-1 

00 

Pd(n) - P^(n) - PnCn-l) - / fo(X) 1 tl-g(X-Xi)] g(X-Xn) ^ (63) 

-00 !■! 
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The mean time to acquire is, 

CO 

T - AT I i Pd(l) 

1-1 

If we choose the set of pointing locations (X^, 1-1 to n) to maximize Pjj(n) , we 
would have n simultaneous equations to satisfy (for X - 1 to n) , 
n 

J f^(X) n |l-g(X-X^)] g^ (x-^) dX- 0 (6S) 

-*■ ■ ■■ 

if‘k 

The values which satisfy these equations could then be ordered to minimize T^ 
given a detection occurred in the n trials. 

These are only two of thn possible cost functions. Whatever cost function we 
choose, one additional factor should be considered* We choose a cost function 
such that a sequence X^, i-1 to n is optimal for n tHals* Then, if the target 
is not detected, we may choose to extend the search for another m trials. Then, 
if the sequence X^, i-1 to n previously defined is contained in the sequence 
X^, 1-1 to n + m which is optimal for n -f m trials, the sequence is uniformly 
optimal. 
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One such uniformly optimal sequence is obtained by simply choosing to maximize 
given the n-1 previous trials. The sequence is obtained by choosing 
then X 2 » etc^ The sequence so defined is not necessarily optimun in terms of 
minimum mem. time to detect the target or maximizing the probability of detection 
in n trials, but it is uniformly optimal. 

The optimal sequence is difficult to define, since we do not know g(X). This 
results from the fact that g(X), the detection probability, is significantly 
affected by atmospheric scintillation. This affect can be modeled by assuming 
the variance of the log of the detected signal (C) is constant, but unknown, and 
that h(G) is the probability density of C. Thus, for our optimal sequence, Xn 
would be chosen to mas^miZe, 


Pd<») - ^ ^ f tl^g(X-^lG)] g(*^lG) h(c> dX dC (66) 

o-«» ° 

1-1 

One more problem remains to be considered. Consider the possibility that a 
cloud was obscuring the target for exactly the first K trials, and let q(K) be 
the probability of this event* Then, 

B-l . . 

,.(n> * E q(k) ' * B [l-g{X-X^lC> B(C) dX « (67) 

K-0 irtthl 

Equation (67) must be solved for the ^ that maximizes p^(n) for n-1 to some 
relatively large number (perhaps up to 200), in order to obtain an optimal 
acquisition search sequence. Then, since the models for q(K) and h(C) are not 
well known, it would probably be desirable to alter the parameters and generate 
a new optimal sequence, and compare the sequences thus generated. Finally, 
g(XlC) is xu)t stattonary with tinw since the range is decreasing rather rapidly 
when acquisition begins. The rate of decrease of range, and hence rate of 
increase in mean signal level, is a function of the distance of the target 
from the Shuttle orbit plane, thus to be as precise as possible, a new optimal 
sequence should be generated for each target pass. 

One minor problem, of course, is that g(XlC) cannot be described in closed form 
when log-normal scintillation is considered. Thus, there is simply no way to 


2-64 


SYNTHESIS AND ANALYSIS OF 

PRECISE SPACEBORNE LASER RANGING SYSTEMS 


REPORT MDC ET729 


FEBURARY 1978 


VOLUME II 


evaluate Equation 67 except numerically, and an enotnxjus computational load is 
anticipated even to generate a abort optimal aequence, 

A more practical approach needa to be deviaed in order to obtain aome baaic 
design data. First* consider the case where the scintillation variance la 
small and the mean signal large. The detection probability will then be 
nearly 1 In the central portion of the beam* and nearly zero outside of the 
central region. Thus* we could approximate. 


g(X) - 1. /X/ < p 
- 0, /X/ > p 


If we arrange a minimum length search sequence so that all possible target 
locations within a specified search radius are within at least one detection 
radius (p), the resulting set of search locations* X^* would be descrlbable 
as a two dimensional sampling grid with any three adjacent points at the comers 
of an equilaterlal triangle of dlmettSlon py|3 on a side* The Overall detection 
probability would be simply the probability that the target was within the region 
searched* assuming to obscuration. This sampling grid structure was chosm for 
the balance of the acquisition studies; the sample spacing was considered a 
variable. 


The mean return signal la a function of the planar distance* v* between the 
beam center and the target* and the e ^ beam radius* r^. 




e 




where a ■ constant detennlned by the link conditions. 


( 68 ) 


When log-normal scintillation is postulated, the observed signal energy, S, is 
a random variable with a probability distribution, 


f(a) 
where C 


1 1 s ^ 

- exp {-[f In C] /2G} 

■ variance of the log of S 


(69) 


Next, assume that the receiver employs an ideal phoco-electron counting, sliding 
window integrator. Then, when the signal is wholly within the window, the 
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probability of detection given Sf for a threshold, L, is, 

K-d 


where S Is the signal energy in photona/pulse 
and n * quantum efficiency of the detector 


Thus, the overall probability of detection Is simply, 

* L-1 . 

T . s. -vT 

K-O e f<S)dS 


1 - 


(70) 


(71) 


This integral was CQnq>uted n\imerically for the limiting case, C • 0.6, for 
tj^ical threshold levels, yielding the curves in Figure 34. The average 
signal is simply nS^, as defined previoualy. An emperical expression was 
derived for these curves, of the form, 

7 r 

1 " i'" ^ 

' K«0 

where the a^ were congiuted numerically for each threshold settings 






(72) 


Next, let the points in the sampling grid be defined as ^ij* 


Xy - id + |- (j mod 2) 




(73) 

(74) 


and 


2 2 2 

hi 

where R la the desired search radius. 


Next, assume that the target is located at a point, }^, Y^. Now, for any 
particular pulse pointed at the 1,^ sample grid location, the distance between 
the beam center and the target is defined by. 
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PROBABILITY OF DETECTION VS AVERAGE SIGNAL 
(WORST CASE ATMOSPHERICS) 




THRESHOLD “ 8 PE/PULSE 


THRESHOLD = 10 PE/PULSE 


THRESHOLD = 12 PE/PULSE 


THRESHOLD * 14 PE/PULSE 







SYNTHESIS AND ANALYSIS OF 

PRECISE SPAGEBORNE LASER RANGING SYSTEMS 


REPORT MDC El 729 

VOLUME II 

FEBURARY 1978 


The average Bigaal expected is, 


■ ri 


w ly! 


<75) 

(76) 


The probability of de'tectioa Isi one coopleta pass through the sampling grid 
ia simply, 

f 10 

P - 1 - E exp ( - E 

iij V K-0 ^ 

The beam radius (r^) has a significant effect on the detection probability. 

If the radius is either very large or very small, the detection probability is 
small, thus we conclude that an optimum beam radius can be defined for any 
specific set of conditions. Figure 35 shows the optimum beam radius as a 
function of the number of samples in the grid and the search region radius, 
where the nominal received signal, from a target on beam center > was 100 p-e/pulse, 
with a 0.5 mrad transmit beam width. 



The detection threshold affects the optimimi value of r^, as shown in Figure 36 
resulting in a decrease of the optimum beam radius as the threshold is increased . 

These calculations were perfo^ -ed under the presumption that the target was 
near the center of the search region. Figure 37 shows the effect of moving 
the target various distances from the center of the search region on the 
probability of missed detection. The effect of varying the threshold on the 
probability of missed detection is shown in Figture 38 . For both of these 
figures, the target was located on a radial inclined 30* from the horizontal 
axis of the saagtllng grid. We found that no significant changes resulted 
from varying the angle of the radial. 


The next step was to assume that the probability distribution of the target 
location was a zero mean circular normal distribution. Thus, the probability 
that the target is at a distance r^ from the center of the sampling grid is. 


p(rj) - Mtp 


( 78 ) 


where. 


a *■ a ■ a • standard deviation of the location error 

2 - 68 . 
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A conservative value of o * 0,2 mrad was assumed, and the overall probability 
of missed detection was computed as a function of the number of samples in 
the grid, as shown in Figure 39. 


WEIGHTED MISS DETECTION PROBABILITY VS NO! OF SAMPLES IN GRID 



Also shown in Figure 39 are three phantom lines, drawn through the origin 
and the lowest points on the computed curves. These lines represent constant 
payoff lines, i.e. , the log of probability of missed detection in linearly 
propoxtlonal to the nund>dr of trials The significance is that the probability 
of missed detection in N trials is the same whether the nimber of trials rep- 
resents maximally Interlaced samples in the grid or repeated trials of a 
smaller slaed grid. Thus* for the most robust link shown, there is no practical 
reason to define a sMipling grid larger than 35 samples. However, for less 
robust links, larger sampling grids represent additional improvement in detec- 
tion probability. The slppea of the bounding lines appear to be linearly pro-: 
portional to nominal signal strength. Thus, we can approximate, for this 
specific set of conditions, 
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(79) 

where N is the number of puleee and E la the expected number of photo-electrons 
per pulse from a target oh beam center where the beam divergence la 0.5 mrad. 

The validity of this approximate e3q>ression can be shown In the following manner. 
GOnsider the case where the target Is located In the center of the sampling 
grid. The sampling grid can be approximated as a set of sampling points located 
on concentric circles about the target. The spacing between adjacent circles 
Is hr. The number of Sample points on the 1th circle from the center Is 
epproxlmately 61. Assume there are n circles, extending to R <> nhr In radial 
distance, with a total number of aamplea, N ■ 1 -F 3n(n-(-l). The mean signal 
strength for a sample on the 1th circle Is, 

*0,1 ’TT* <8Q) 

■ o 

As noted previously, the probability of detection for this sample Is, 

/“ A 

P* 1 - exp - 1 2- . A,, S, j J 




Thus, the log of the probability of missed detection for the 61 trials on the 
1 th circle Is, 


-61 Eajc/S') "o( 

K-1 Vo / 


-K(lAr/r^j 
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(82) 


(83) 


Thus» to « first order approxisiatloef the log of the probability of missed 
detection Is lines. :ly proportional to both K and signal energy per pulse » 

and Inversely proportional to the acquisition search area. As long as the 
search radius is considerably larger than the standard deviation of the target 


location uncertainty, and the target grid spacing Is small compared to the 
nominal beam radius, this approximation should be quite reasonable regardless 
of target location within the grid, as noted during the numerical analysis. 
Furthermore, note that r^, the beam radius for acquisition, does not appear 
^ a first order effect on che mlsaed detection probabll Lty, provided that 
r^ Is large enough that the second and higher order terms do not effect che 
results. 


Thus, we conclude that for reasonable system parameters, acquisition is likely 

to be quite tellable even for worst case atmospheric scintillation. Further, 

note that at nominal acquisition conditions , i.e>, 20^ elevation, the range 

is reducing at 'V' 7.5kth/sec, thus the mean return signal energy will lncrea!ie 

almost 40% in Che first 10 seconds of the acquisition process, further enhancing 

the acquisition probability. Also, we note that the Optimum beam radius for 

acquisition Is sufficiently large that che second and higher ordered terms 

in the approximation are negllgable. Thus, we find that choosing a beam radius 

somewhat larger than the calculated optimum yields better results than choosing 

a smaller beam radius (the are alternating In sign) . Finally, we found 

that the controlllhg paraineCers were the nominal energy per pulse and the grid 
2 

sample density (R /N) . 
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2.4 Stable Platforms 

A number of experiment pointing systems have been proposed for Shuttle/Spacelab 
applications. The general thrust of these pointing systems is to provide a 
general purpose mount which is capable of pointing experiments with varying 
levels of accuracy and stability. Four of these systems were evaluated for 
compatibility with a laaer ranging experiment. The initial screening resulted 
in concluding that only one of the proposed pointing systems was potentially 
suitable for a laser ranging experiment. The four pointing systems considered 
were the Instrument Pointing System (IPS), the HLuiature Pointing Mount (MPM) , 
the Small Instrument Pointing System (SIPS), and the Annular Suspension and 
Pointing System (ASPS) . 

Both the IPS and the MPM employ mechanical isolators to reduce the sensitivity 
of the pointing system to Shuttle disturbances. These Isolators allow the 
glmbal centroid to move with respect to the Shuttle structure. Any unknown 
motion of the glmbal centroid with respect to the Shuttle coordinate system 
results In a residual translation error of the ranging data. Consequently, 
both the IPS and the MPM were considered basically unsuitable for a precision 
laser ranging experiment. 

The ASPS concept includes a magnetically suspended fine pointing glmbal system 
Inside a conventional coarse glmbal system. This concept was not considered 
suitable for applications demanding significant angular agility. 

The remaining pointing system, SIPS, Is considered basically suitable for a 
laser ranging experiment; no fundamental incompatibilities were identified, 
although the design slew rates and angular acceleration capabilities will 
limit the experiment capability to make ranging measurements to target grids 
composed of many and/or widely spaced targets. A number of areas of concern 
were Identified, as discussed In the following paragraphs. 

The first area of concern is one of safety* This concern is due to the design 
of the SIPS gimbals, which allow full hemispherical pointing capability. 

Since the laser ranging system employs a high peak power laser transmitter, 
a system of interlocks and control procedures will be required to ensure that 
no significant amount of laser energy is incident on any part of the Shuttle, 
and spectrally selective window covers are recommended for all windows in the 
Shuttle and Spacelab which could be, or could view any part of tne Shuttle 
illuminated by the laser ranging system transmitter. 
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The second area of concern reaolCB frou the requlrentent to transfer the range 
measurement to the Shuttle center of gravity (eg). This transfer requires 
knowledge of the SIPS glmbal angles and position of the glmbal centroids In 
Shuttle coordinates. As discussed In Section 2.2, the location of Che Shuttle 
eg Is Included in the u.atheaiatlcal model of the system. However, any m\certaluty 
In glmbal angles or glmbal position resitlts in a coupled ranging error. As 
currently Implemented, the SIPS glmbal tmgles are sensed with 12 bit precision 
shaft encoders (vl,5 mrad resolution). 

The distance from Che glmbal centroid to Che Shuttle eg Is a function of the 
location of the experiment within the payload bay, and can range from 2.5 to 
12 meters; the maximum ranging error coupled Into the measurement by glmbal 
angle readout resolution Is on the order of 1.5 x 10 ^ x h/^2 (rms). 

Note Chat angular deflections due the thermal gradients and/or structural 
flexures could add directly to this imcertalnty, but for the purposes of 
this study were assumed to be an order of magnitude smaller than the readout 
resolution. 

The third iirea of concern is the relationship between Shuttle attitude, during 
the pass, and the required angular accelerations about the glmbal axes to keep 
Che pointing vector oriented coward the target. This la a concern b«=Gau8e of 
the limited ^mgular travel of both the inner (right-left) glmbal aid the eleva- 
tion glmbal. The Shuttle orientation strategy devised in the mission simulation 
(Section 2.2) resulted in si'lov*tiug a modified, x-axls perpendicular to orbit 
plane, I POP) liuutial attitodo such th.U the Shuttle Z axis intercepted 
the target grid center at the point of closest approach. The target grid, 
under this condition, wiis located within a few degrees of the Shuttle Y-2 
plane throughout the ranging pas^i. Tills is not mi optimum attitude strategy 
for the SIPS mounted laser ranging experiment, since either an aalmuch glmbal 
flip is required, or the SIPS elevation axis design would have to be modified 
to permit a 180^^ elevation capability. The simplest solution seems to be to 
tilt Che Shuttle Z axis by an angle, 1, from the grid center at the point of 

closest approach (|•etalnlng the Y .axis in the orbit plane). Then, the maximum 

2 

angular av^oieratiou required about the azimuth axis is approximately . 02/sin*' 

(degrees/sec*') , Cor a Shuttle altitude of v330 km (180 uml) . The maximum azimuth 

2 

glmbal acceleration capability is assumed to be 0.2" /second”, thus a tilt angle, 
1, of at least 18.5" Is needed Co allow continuous tracking at worst case condi- 
tions, It is somewhat undesirable to choose considerably larger tilt angles. 
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since this would bring the laser beam closer to the Shuttle or Spacelab structure, 
which is undesirable from a safety standpoint. 

The fourth area of concern Is the Interface between the laser ranging experi- 
ment and the SIPS computer. As currently defined, the experiment will derive 
permitting error signals In up-down, and left-right coordinates which are 
forwarded to the SIPS computer to Initiate corrective action. The concern is 
related to the effect of lags caused by data transfer delays and finite computer 
cycle times on the pointing error. Consider the following scenario. The 
experiment computer performs a number of essential functions, including a naviga- 
tion function and an attitude reference fwction. The navigation function Is 
Included so that ranging data, after initial target acquisition, can be used 
to update the spacecraft position data and reduce the net pointing uncertainty 
for the balance of the ranging pass. The navigation data Is used to predict 
the proper pointing direction and to position the range gate for the next scheduled 
transmit pulse. We anticipate this prediction will include the target range 
vector and the range rate vector. In stellar Inertial coordinates. The attitude 
reference function maintains a current estimate of the actual pointing direction 
of the laser ranging experiment. The pointing command function is derived from 
the current pointing direction and the desired pointing direction for the next 
transmit pulse time. A prediction algorithm generates a series of pointing 
commands which will result in being at the correct pointing direction at the time 
of the next scheduled transmit pulse, with the correct angular velocity to main- 
tain track during the pulse propagation time. This time period varies from two 
to six milliseconds, and the net change of direction during this Interval is 
2 V^/C, also referred to as the point ahead angle. At the point of closest 
approach, the point-ahead angle is about 55 prad, thus the maximum rotation 
rate of the pointing vector is lees tha .1 10 prad/ms. If the SIPS Gonqputer can- 
not accept and process a pointing error command in a prompt and precisely r ce- 
dlctable manner, a dynamic error is Introduced, which is proportional to 
the variation of the execution time and the pointing vector rotation rate. 

The areas of concern discussed above are not insuperable problems; we conclude 
that a laser ranging experiment can be configured using the SIPS to point the 
transmit and receive telescopes. Integration of the laser ranging system and 
the SIPS, however, is considerably more complex than envisioned for the normal 
Instrument payload for a SIPS application. 
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3oO Experiment Preliminary Design 


3.1 Pallet Mounted Experiment 

Hie configuration of the pallet mounted experiment is Illustrated in Figure <i0. 
The experiment hArdware is composed of two groups, the support electronics, 
mounted on coldplates attached to the pallet side, and the optical bench and 
associ-iLed components, which are mounted in a protective enclosure on a sill 
level platform. A removable cover is provided to shield the optical components 
from contamination and to maintain thermal stability when the experiment is 
not operating. 

The optical bench mounted components include the laser transmitter, the optical 
receiver, star trackers, the attitude reference gyros, and the telescope and 
beam steering mirror. 

The active elements of the system which are mounted on the optical bench are 
thermally Insulated from the bench; energy which Is dissipated within these 
components is discharged through an active thermoelectric system to the pallet 
supplied liquid loop. The enclosure is thermally insulated and equipped with 
heaters to maintain a stable thermal environment. 

The support electronics components are mounted on coldplates connected to the 
pallet supplied liquid cooling system. If required, a thermoelectric mounting 
plate can be inserted between the component baseplate and the coldplate, 
allowing a precise bidirectional heat flow control to maintain a narrow 
temperature band within the component. 


Figure 41 pictorlally describes a typicei ranging experiment. Approximately 
45 ntlnutes prior to the beginning of the rangitig pass, the Shuttle is maneuvered 
to the nominal Inertial orientation for the pass, and the experiment attitude 
reference system is initialized and calibrated during the next 15 to 20 minutes. 
The Shuttle orientation .s nominally X-axis perpendicular to orbit plane (X-l’Or) , 
with the Z-axls pointed directly at the center of the target grid at the point 
of closest approach. The Shuttle angular rates are controlled to as low as 
possible to minimize the need for thruster activity during the actual ranging 
operation. Wlien the first target comes into view (^ 120 ” elevation at the ground 
target), an acquisition sequence Is initiated. A short target verify period 
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follows c‘ie Initial target deteetion. The experiment tlien eaters the nominal 
ranging mode, whieh continues until the last target is no longer in view. At 
this time, the experiment is shut down and normal Shuttle activities resume. 


LASER RANGING EXPERIMENT PICTORIAL 

ACTIVE 



ALIGNMENT 

Figure 41 

The data accumulated during the ranging pass is recorded for subsequent trans- 
mission to the ground. A typical individual measurement record includes both 
the transmit pulse departure time and the time of arrival of the reflected 
signal, attitude and giinbal angle data, a target label, the current Shuttle 
position data, and miscellaneous system parameters pertinent to data reduction. 

3.1.1 Functional Description 

The major elements of the laser ranging system on’-board the Shuttle are shown 
in Figure 42. The major interfaces are aVmwn to illustrate the move signi- 
ficant functions. The clock subsystem issues a command sequence to the laser 
transmitter to initiate a pulse at cr near an expected time, nominally at a 
10 pps rate. The clock also determines both the actual pulse dep.irture time 
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and the detected return pulse time, and forwards this data to the data assem- 
bler and the navigation computer. The navigation computer either accepts or 
rejects the data point based on several test criteria. If accepted, the 
measurement is used to update the state vector via a Kalman filter technique. 

The navigation computer also generates the spaeecraf t/target range vector (and 
its derivative) for the next pulse emission time. The magnitude of the range 
vector la used to set the ranging receiver range gate. The range vector and 
its time derivative is forwarded to the ARS computer, which rotates both into 
Shuttle coordinates and converts this data into gimbal position eotrattand. The 
navigation computer outputs range vector data at a 10 sps rate; the ARS computer 
processes the direction and rate predictions into a gimbal angle trajectory 
command sequence which is delivered to the pointing control subsystem computer 
at a nominal 200 sps rate. Tlie pointing control system compares the gimbal posi- 
tions with the command positions and issues the necessary acceleration commands 
to the beam steering subsystem torquers. The ARS computer also maintains a 
current estimate of the optical bench attitude In stellar Inertial coordinates. 

The laser ranging system Is controlled by the cotimand and control system which 
contains the desired ranging measurement sequence and the necessary navigation 
data to initialize the navigation computer prior to initial target acquisition. 

The RAU is the data and command interface with the Shuttle/Spacelab data and 
command and control system. Shuttle electrical power is conditioned and 
delivered to the laser ranging system components by the power conditioning unit. 

The thermal control subsystem includes Insulation, heaters, coldplates, and 
thermoelectric heat transfer devices to maintain the components within the 
desired operating temperature ranges. 

Table 4 summarizes a typical laser ranging exi^criment operation during a 
ranging pass. The first operation, called pre-allgn, requires applying power 
to some of the active components, including the thermal control subsystem, to 
establish thermal equilibrium in the sensitive elements. Tlie next step is to 
Initialize the navigation computer and commence the navigation function. This 
includes determination of the desired Shuttle attitude for the ranging pass. 

Tlie attitude data (stellar Inertial coordinates) is forwarded to the Shuttle, 
which is rotated into the desired attitude for the pass. 
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Table 4 


TYPICAL LASER RANGING EXPERIMENT 


Mode 

Activity 

Duration 

Control 

Prealign 

Powerup Selected 
Components 

TBD 

External Command 

Navigation 

Initialize 

Nav* Computer on 

TBD 

(1) External Cmnd On 

(2) Transfer Desired 
Shuttle Orientation 

ARS Align 

ARS on 

Star Trackers on 

'^'30 min 

Automatic 

Coast 


'Vis min 

Automatic 

Pre-Acquisition 

Gimbals Active 
Laser On 

*^20 see 

Automatic 

Acquisition 

Dither about Nominal 
Target Position 

'vlO see 

Preprogrammed 

Sequence 

Verify 

Test & Update 

*v2 sec 

Automatic 

Ranging 

Track/ Slew/Track 

'vZbO sec 

Preprogrammed 

Sequence 

Post-Pass 

Laser Off, Thermal 

'V'S min 

Automatic 


Cooldown 

The next step Is to enter the ARS align mode. The star tracker output data 
Is used to determine the actual Inertial orientation of the optical bench and 
to calibrate the ARS gyros drift and scale factors. This mode Is maintained 
until the star trackers line of sight is (effectively) obscured by the earth. 

At this time the ARS reverts to a gyro only mode until the end of the ranging 
pass . 

The next phase Is the preacquisition mode. In this mode, thei navigation 
computer begins the target view computation sequence, and the ARS computer 
uncages and commands the desired pointing angles for the mirror gimbals. 
Approximately 20 seconds before the start of the acquisition sequence, the 
laser transmitter is enabled and thermally stabilized. During this period, the 
mirror is positioned to the nominal predicted location of the first target 
to be observed. 
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IC the target Is not detected during the preacquisition mode, the system 
enters the acquisition mode, which requires dithering the mirror gimbals about 
the nominal line of sight to the target In a preprogrammed sequence > 

Once the target Is detected, the system enters the target verify mode. In this 
mode, the observed detection characteristics are compared to preset criteria, 
and a valid/invalid decision is made. If a valid target decision is made, the 
navigation computer is updated and the system enters the ranging mode. If, 
however, the invalid target decision is made, the acquisition mode is resumed. 

An alternate acquisition mode was considered which merged the acquisition and 
target verify modes. In this alternate eoncept, a sequential search was per- 
formed In space about the nominal line of sight. 

However, Instead of moving from point to point In the search sequence after a 
prescribed dwell time, each point would be viewed until either a no target or 
a target present decision could be made with a p^rescribed confidence level. 

This acquisition mode was not considered as suitable for this appllGation as 
the fixed dwell time concept, since the major acquisition problem is due to 
atmospheric scintillation, and deep fades can last a considerable time period. 

In the ranging mode, a preprogranimed measurement sequence is followed. The 
transmitter is slewed to illuminate one of the targets in the grid, and ranging 
measurements are made for N successive pulses, where S is a variable dependent 
on the number of targets in the grid. 

The system is then slewed to the next target In the sequence and the measure- 
ment sequence repeated. This sequence of track, slew, track is continued 
until the last target is no longer in useful view. 

The post-pass mode is provided to ensure completion of all buffered functions 
and CO ensure that the laser has had time to cool dot«m. The laser ranging 
system is subsequently powered down to the hold mode, where only the clock and 
TCS monitor functions are retained. 

3,1,2 Folding Mirror and Pointing Control System 

The baseline model for the folding mirror, gimbals, and control system is 
derived from a similar system under development for an Air Force laser coiranuni- 
catlon program. Figure 43 shows the overall dimensions of the glmbal mirror 
assembly. It is in a 2 axis roll-elevation configuration with elevation inside 
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The gimbals are eleecromagnetically caged (during launch) with power off, 
and can be uncaged by application of power to redundant D.C. torque motors. 

Figure 44 shows the assembly details. 

The gtmbal control system is digitally implemented, excluding the torque motor 
power amplifiers* Figure 45 shows one axis of the glmbal control system. 

The normal linear control system is augmented with an acceleration schedule 
control to facilitate large angular steps. The acceleration schedule control 
applies maximum acceleration torque in the desired direction until the error 
signal is less than one half of the initial commanded change in pointing angle. 

At this time the commanded acceleration torque is reversed The acceleration 
schedule Is terminated when either the error signal changes sign or a computed 
time has elapsed. The (conceptual) switches are then positioned to enable the 
normal linear mode control function. 

The control system from which this model is drawn has been demonstrated in 
both simulation and brass-board testing to attain '^2 yrad rms open loop 
pointing stability but the design application required only very slow 
slewing rates. The major change to the reference model, therefore, is to adapt 
to our requirement for minimum time to slew from one target in a grid to another. 
A number of factors must be considered in order to establish the feasibility 
of these changes and the expected time to settle after a sizable angular step. 

The first concern is the ability of the inductosyns to follow relatively large 
angular rates. The inductosyn is an incremental counting device which maintains 
the 20 bit accuracy by counting changes (+) in the least significant bit. The 
subsystem presently used has an error free maximum angular rate of ^^0.06 rad/s 
(3,4*/s). If the angular rate exceeds the error free rate, the output 
lags the shaft angle position. The maximum lag is ^^700 urad. If the angular 
rate exceeds about 7 rad/s (400“/s), a loss of reference occurs and the 
gimbals must be stopped to reestablish the actual position. 

2 

The maximum angular acceleration for the existing system gimbals is rad^s“ 

2 

in elevation and 2 rad/s in roll. The maximum angular step for the elevation 
glmbal is on the order of 0.7 rad, thus the peak angular velocity would be on 
the order of 2.4 rad/s, well below the angular velocity at which loss of 
reference would occur. The maximum angular step in the roll axis is less than 

* MDC Report E1492, Page 127 
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Inductosyn, the angular position is a direct readout, thus when the angular 
velocity decreases below the 6“/s rate, full angular position aecuraey is 
immediately restored). At angular rates above b“/s, the optical shaft 
encoder output becomes increasingly inaccurate, with the error doubling for 
each factor of two in angular rate. 

Tims, we conclude c!»at the most significant problem with either type of 
encoder is that our ability to determine when the gimbal has reached the 
half-way point is decreased, compared to static condtions, at large angular 
velocities. 

The next question to be addressed is the control technique to be employed. 

The reference system uses an all digital control system Implemented with a 
microcomputer. Consequently, a change in control law simply requires 
programming a new memory chip. 

The reference control system employs a 5 ms cycle time. The cycle time can 
be reduced, if necessary, at the expense of increased power consumption. 

At a maximum angular velocity of 'v2.4 rad/s, the position of the gimbal would 
change v,12 mrad between successive samples, tluis we see that the loss of 
accuracy due to output lag for the indue tosyn (or resolution for the optical 
encoder) is trivial in comparison. 

As presently planned, tne acceleration schedule changes sign when the error 
sigtuil is less than one-half of the commanded angle change. Thus, there is 
a lag between the actual hal f-anvle rrns.sino and Hiia r'nmmpTirement of the 
negative acceleration schedule of up to one cycle time. 

This lag results in a residual angular velocity at the end of the acceleration 
schedule, virtually Insuring an overshoot condition. The residual angular 
velocity is a function of both the commanded step size and the cycle time. 

This residual angular velocity increases the time required for the gimbal to 
settle to within acceptable pointing accuracy (on the order of 50 ^trad) . 

The settling time can be quite long if the residual velocity causes the 
gimbal to significantly overshoot the linear control range. 

The next concern Is to examine the situation where the target gimbal angle 
Is not stationary. Consider the situation where the Shuttle has a non-zero 
angular rate In the gimbal axis coordinates. Thus, the commanded gimbal 
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angle ie a time varying angle defined as 0^ “ An Ideal acceleration 

schedule vould accelerate to time and then decelerate till time the 
glmbal angle ai\d angle rate wovild match the target angle and angle rate if. 



where the sign of n 1 b' chosen* to make t^ (hence t,>) positive and real. 

When the angular step slat' is large and the target angular rate modest, we 
find that the error sigival time Is, 

E " lW2 + 0 “ /4a 
o c 

Consequently, we see that the half-angle step strategy results in jj small 
switching lag (cowv>ared to an ideal acceleration schedule). 

There is aiiother factor, however, which adds to the xitility of the half ai\gie 
concept. The torque actually a\ipiled to the gimbals for a fixed drive voltage 
varies somewhat with glmbal angle. This torque ripple results In a potentially 
significant dispersion of angular position and velocity at the end of a time 
ordered acceleriition profile*. The half angle swltchlvig strategy lends to 
minimise this dispers'on. However, it is this dispersion which requires 
adding the timer to terminate the acceleration schedule, since a greater than 
expected torqiu-* during the latter half of the schoviule could result in can- 
cellation of the lag effects to the extent that no overshoot occurs. 

A simplified simulation was cmisCructed to test the perfomauce of the glmbal 
control system. Torque ripple was modelled as a random process; encoder error 
was also modeled as a random process with an angular velocity dependent 
variance* Figure 46 shows the response ot the roll glmbal to a 0.13 rad step 
cptmnand with aero ai\d + 50 mrad/s command rates. The glmbal settled to 
within 50 prad within approxinpitely 0,7 s. The cvnitrol loop was I'peraled 
with a S ms cycle time. Figure 47 shows similar curves for the inuev glmbal, 
also at a 5 ms cycle time. In this figure, more pronounced overshoot can 
be observed; oven though the first crossing occurred much sooner than in 
the roll glmbal case, the overall time to settle was only slightly less 
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«5S s)* Figure A8 shows the response of the inner glmbal to a 150 mrad 
step with + or - 50 mrad/s rate, but with a 1 ms cycle tlmep The large 
overshoot observed In the preceding figure has been eliminated* and the gimbal 
settled to within 50 mrad In 'v 0.35 s. 

In this simulation, the gimbals were modeled as frictionless and linear. 

This assumption does not materially affect the dynamic performance except 
in the small angular rate regime where breakout forces are encountered. 

In conclusion* a modified version of a pointing control system currently 
under development has been evaluated for the laser ranging experiment, and 
found to be reasonably close to meeting all pointing requirements. Add!*- 
tional modifications to increase the avalxable torque and decrease the time 
to slew and settle are certainly plausible; considerable additional evalua- 
tion would be required, however* to verify the fine pointing capability if 
larger torques are used. 

3.1.3 Opto-Me chan leal Subsystem 

The opto-mechanlcal subsystem Is shown schematically in Figure 49. The 
laser output Is passed through a frequency doubler and then expanded to approxi- 
mately 2 cm diameter and defoeussed slightly to obtain 0.5 mrad divergence. The 
doubled output 1s reflected by a dichrole mirror, and directed to an alignment 
mirror, which directs the beam to the final folding mirror. The final folding 
mirror folds the transmit beam to be aligned with the primary telescope axis . 

A small portion of the transmit beam Is extracted at the final folding mirror 
and used to determine the overall transmit beam alignment. 

Another small fraction of the transmit energy is diverted through a neutral 
density filter and reflected to the high speed photomultiplier tube in order 
to obtain an estimate of pulse departure time. 

The receive optics path Includes a bifurcating mirror which allows a star 
sensor to be used to calibrate the gimballed folding mirror (not shown) . The 
central part of the receiver field of view passes through the bifurcating mirror, 
encounters two torque motor beam steering elements, and a narrow band optical 
filter. Provisions were made to insert either a shutter or a variable density 
filter. The receive path Includes another bifurcating mirror, which splits off 


3^16 


GIMBAL ANGLE (MRADJ 


RESPONSE OF INNER GIHBAL 



TIME (SEG) 













SYNTHESIS AND ANALYSIS OF 

PRECISE SPACEBORNE LASER RANGING SYSTEMS 


QMS SCHEMATIC 


REPORT MDC El 729 
FEBURARY 1978 


VOLUME II 



the outer edges of the Inner field of view to a quadrant PMT for gross receiver 

boresight error coinpensatlon. The central 0.5 mrad field of view is then 

focused on the ranging detector. 

The collecting and resolving power of the telescope is such that the star sensor 

path could be split into two paths (50X reflective mirror), with the second 

path directed to a television camera to enhance the experiment monitoring 
capability of the crew. 

Various concepts for on-orbit calibration were explored to determine relative 
accuracy and feasibility. The first phase of the on-orbit calibration process 
would be to calibrate the glmballed folding mirror and receiver primary focal 
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plane using the focal plane star sensor. The attitude reference system must be 
operating during this period to relate star sensor output to the navigation 
base attitude. The boresight location of the inner field of view can be deter- 
mined by driving the gimbals such that a star traverses through the inner field 
of view in a prescribed manner, thus the star sensor drop out during the 
traverse is indicative of the centering of the inner field of view and the star 
sensor field of view. 

If the narrow band optical filter is removable, and the PPMT is replaced with 
a CCD type star sensor, a further fine resolution calibration of the torque 
motor beam steering units can be performed, resulting in a very precise 
centering of the inner field of view with the outer field of view and the 
gimbals. 

The major problem remaining is to verify the alignment of the transmit beam 
with the calibrated receiver boresight. The most attractive approach seems to 
be to locate a collimator on the optics bench, and to command the ginfcal system 
to a preselected position such that the transmit beam illuminates the collimator. 
Of course, this would be beyond the region where the glmbal calibration can be 
verified by stellar observation, but the probable error resulting from eKtrapola- 
tion is likely to be much smaller than the alignment tolerance. Once this 
calibration is completed, the alignment sensor is also calibrated, and subse- 
quent short terra alignment variations can be removed using the alignment mirror 
adjustment capability. 

The receiver and gimbal calibration can be readily accomplished during the gyro 
drift and scale factor calibatlon period preceding a ranging pass. The trans- 
mitter calibration can be performied in the short time inteirval between initia- 
tion of laser operation and the start of the ranging pass. The transmitter Is 
turned on and allowed to stabilize, and then the alignment verified (or 
corrected). The gimbals are then commanded to the nominal acquisition atti- 
tude and normal operation begun. The location of the collimator is chosen to 
minimize the gimbal slew required to move to the acquisition orientation. 

3.1,4 Interfaces 

Tlie experiment interfaces are with the pallet and pallet subsystems, and 
ultimately with the Spacelab and Shuttle subsystems. The optical bench is 
supported on a sill level platform to obtain an adequate field of view for 
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the ranging experiment. The balance of the laser ranging system components 
are mounted on coldplates attached to the pallet structure, as sho\m in 
Figure 40. 

Most of the LRS components can be mounted directly on coldplates serviced by 
the pallet freon loop. A few components, such as the laser transmitter, will 
require more precise thermal control than provided by the standard pallet 
thermal control system. For these components, a thermoelectric unit is inserted 
between the coldplate and the unit. Heaters and Insulation are provided to 
maintain the optical bench at a stable temperature. 

The LRS data interface with Shuttle/Spaeelab is through a standard As 

presently plivnned, data will be transferred to the Shuttle/Spaeelab in several 
manners. Ranging data will be transferred in discrete packages, assembled by 

the LRS. Each package will be on the order of two hundred bits, and will be 

output at a ten packages per second rate. Performance analysis data will also 
be forwarded in discrete packages, on the order of 5000 bits per package, also 
at the ten package per second ra.te. Engineering data will be assembled and 
forwarded at a nominal 400 bps rate during the operating period, llie data load 

is summarized in Table 5, Ranging data is the data necessary to reduce the 

range measurements on the ground after the pass. Performance analysis data is 
intended to give as complete a picture of the actual experiment performance as 
possible, and may be monitored in part by the experiment operator. Engineering 
data is Intended to monitor the performance of the support subsystems, and is 
expected to be monitored by the experiment operator, as well as being forwarded 
to the ground. 

Table 6 summarizes the expected data load from a single LRS ranging pass. 

Table 7 summarizes the data load if constrained to Shuttle standard data rates. 
It may not be necessary to dump the performance analysis data at the LRS pulse 
rate, which could reduce the peak data rate and total data load considerably. 

Another interface with the Shuttle/Spaeelab is the LRS input data and command 
interface. Tills data Is summarized in Table 8. 

Electrical power received from the Shuttle/Spaeelab Is conditioned and delivered 
to the various components. The laser power supply is the single largest load, 
about 390 watts prime power Is the estimated requirement. The thermoelectric 
coolers will require only modest amounts of power, except under extreme opera- 
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TABLE 5 


DATA LOAD SUMMARY 


I RANGING DATA SUMMARY 


BIT RATE 


FUNCTION 

(NO.) 

(BITS) 

ACQ 

PASS 

TARGET LABEL 

1 

8 

80 

80 

ACQ. GRID LABEL 

1 

8 

80 

80 

Tx 

1 

54 

540 

540 

Tr 

1 

54 

0 

540 

SIGNAL STRENGTH 

1 

12 

0 

120 

GATE WIDTH 

1 

16 

L60 

160 

GIMBAL READOUTS 

2 

y 

380 

380 

SUBTOTAL 


/ 

1240 

190G 

OVERHEAD 


/ 

138 

211 

TOTAL 


/ 

/ 

1378 

2111 


II PERFORMANCE ANALYSIS DATA 


BIT RATE 


FUNCTION 

(N^f 

(BITS) 

ACQ 

PASS 

STATE VECTOR 

X 13 

32 

4160 

416G 

• COVARIANCE MATRIX 

X 91 

32 

0 

29120 

RANGE ESTIMATE 

1 

32 

320 

320 

~ ~ "'RANCS'"RES-XPUAIj -a^.w*’***"'*^ 

1 

16 

0 

160 

DRAG ESTIMATE 

1 

16 

160 

160 

ATM CORRECTION 

2 

16 

320 

320 

TARGET POINTING VECTOR 

3 

32 

960 

960 

TARGET POINTING RATE 

3 

36 

480 

480 

QUAD. DETECTORS 

4 

8 

320 

320 

BASEPLATE ORIENTATION 

3 

32 

960 

960 

GIMBAL ACCEL. CMNDS 

2 

8 

3200 

3200 

subtotal 



10880 

40160 

OVERHEAD 



1209 

4462 

TOTAL 



12G89 

44622 


ENGINEERING DATA 


FUNCTION (NO.) 

GMT .1 

TEMPERATURE 100 

VOLTAGE 100 

PRESSURE 10 

DISCRETES 20 

EVENT TIMES 1 

ARS DATA (EXP.) 3 

ARS DATA (SHUTTLE) 3 

INS DATA (SHUTTLE) 6 

THRUSTER FIRING LABEL 1 

THRUSTER FIRING TIME 1 

THRUSTER FIRING DURATION 1 

SPARE CHANNELS 
SUBTOTALS 
OVERHEAD 
TOTAL 


(BITS) 

(SAMPLE RATE) 

BIT PATE 

24 

1/60 

0.4 

8 

0.1 

80 

8 

0.1 

80 

8 

0.1 

8 

1 

1.0 

20 

32 

0.1 

3.2 

24 

0.1 

7.2 

24 

0.1 

7,2 

20 

0.1 

12,0 

16 

0.1 

1.6 

24 

0.1 

2,4 

12 

0.1 

1.2 

136.8 

360 



40 



400 
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TABLE 6 




DATA LOAD 

SUMMARY 



TOTAL 



I 

II 

III 

TOTAL 

BITS 

MODE 

DURATION 

BPS 

BPS 

BPS 

BPS 

MBITS 

PRE-PASS 

45 MIN. 

0 

0 

400 

400 

1.08 

ACQ. 

60 SEC 

1378 

12089 

400 

13467 

0.808 

RANGING 

260 SEC. 

2111 

44622 

400 

46733 

12.15 

POST-PASS 

5 MIN. 

0 

0 

400 

400 

0.12 


TOTAL DATA LOAD/PASS 


14.16 




TABLE 7 




DATA LOAD- 

-STANDARD SHUTTLE BIT 

RATES 


MODE 

DURATION 

ACTUAI, DATA RATE 
BPS 

STANDARD 
BIT RATE 
BPS 

TOTAL DATA 
LOAD 
MBPS 

PRE-PASS 

45 MIN. 

400 

400 

1.08 

ACQ. 

60 SEC. 

13467 

25600 

1.536 

RANGING 

260 SEC 

46733 

51200 

13.312 

POST-PASS 

5 MIN. 

400 

400 

0.12 

TOTAL DATA LOAD/PASS 



16.048 
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TABLE 8 

COMMAND LOAD SUMMARY 


GROUND TO LRS 


ITEM 

BITS/WGRD 

NO. OF WORDS 

NO, BITS 

WAKE UP TIME 

24 

1 

24 

TARGET LABEL 

8 

9 

72 

TARGET LAT. 

32 

9 

288 

TARGET LONG. 

32 

9 

288 

TARGET COV. MATRIX 

32 

45 

1440 

ACQ. SEQUENCE 

8 

200 

1600 

TOTAL 



3712 

SHUTTLE/SPACELAB TO LRS - 

INITIALIZE 



ON /OFF /MODE 

8 

10 

80 

EPHEMERIS 

32 

6 

192 

EPOCH 

24 

1 

24 

ARS 

24 

3 

72 

GMT 

24 

1 

24 




392 

LRS TO SHUTTLE/SPACELAB. 




ATTITUDE COMMAND 

24 

3 

72 

MODE/DATA RATE 

8 

1 

8 


ting conditions (max. coolant temperature) . Assuming the laser is to be main- 
tained at a 20“C heat exchanger input condition » the peak 40“G coolant tempera- 
ture will require pumping '\'300W up 20®C. At that teiiq)erature differential, a 
single stage cooler pumps at about a 1.75 ratio, thus the cooler would require 
about 170 watts to pump 300 watts across a 20°C differential. At the freon 
system lower temperature level, the cooler drive can be reversed to niaintain 
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a steady 20®G heat exchanger input. The overall experiment power consumption 
at full operating condition is estimated to be 854 watts. The estimated energy 
consumption profile is shown in Table 9. About 6.7 kw-hrs are consumed by 
the frequency standard and clock if operated continuously. The energy consumed 
per pass is about 0.42 kw-hrs» most during the 45 minute alignment period 
Preceding the actual ranging pass . For planning purposes , it was assumed that 
ten ranging passes, total, would be exorcised to cover the ground truth sites 
and to calibrate the system, resulting in a net energy consumption of 10.9 kw- 
hours. Table 10 summarizes the estimated weight and power consumption and 
dissipation of the various subsystems and components . The ARS sensors and the 
laser power supply are expected to Interface directly with the Shuttle/Spacelab 
power bus. All other components are powered by the pcx^er conditioner, which 
converts Shuttle/Spacelab bus power into regulated power at the appropriate 
voltages. 


TABLE 9 

POWER PROFILE ESTIMATE 


POWER VS MODE (DURATION) 


ENERGY 

CONSUMED 


PRIME POWER 
INPUTS 

OFF 

(7 DAYS) 

ALIGN 
(45 MIN) 

PASS 

(280 SEC) 

POST-PASS 
(5 MIN) 

10 PAS£ 
(KW-HI 

LASER PWR SUPPLY 

0 

0 

390W 

0 

.3 

ARS SENSORS 

0 

66W 

66W 

66 

0.6 

POWER CONDITIONER 

40W 

359W 

398W 

359W 

10.0 

TOTAL 

40W 

425W 

854W 

425W 

10.9 

ENERGY - 10 PASSES 

6,7 

3.2 

0.7 

0.3 

10.9 
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TABLE iO 

EQUIPMENT CHARAGTERISTICS SUMMARY 


ITEM 

WEIGHT’ ESTIMATE 
(LBS) 

DISSIPATED 

POWER 

j;WATTS} 

COOLING 

POWER 

(WATTS) 

PRIME 

POWER 

(WATTS) 

ARS SENSORS 

34 

6b 

0 

66 

LASER TRANSMITTER 

20 

300 

0-170 

0 

LASER POWER SUPPLY 

10 

90 

0 

390 

FREQ. STANDARD 

20 

20 

0-11 

0 

CLOCK 

10 

10 

0 

0 

OPTICAL RECEIVER 

30 

10 

0 

0 

POINTING CONTROLLER 

25 

20 

0 

0 

GIMBAL MIRROR ASSY 

35 

0 

0 

0 

OPTICAL BENCH 

18 

0 

0 

0 

NAV. COMPUTER 

7 

15 

0 

0 

ARS COMPUTER 

7 

15 

0 

0 

C&C UNIT 

5 

5 

0 

0 

DATA ASSEMBLER 

5 

10 

0 

0 

subtotals 

226 

561 

0-181 

456 

TCS 

15 

20 

G 

0 

POWER CONDITIONER 

25 

92 

0 

398 

TOTAL 

266 

673 

181 

854 
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3.2 SIPS Mounted Experiment 

The configuration of the SIPS mounted experiment Is Illustrated in Figure 50. 

The experiment hardware is shown in two basic groupings, those components 
mounted on the optical bench, and those components attached to the SIPS 
cannister directly. The optical bench is centered in the cannlster, and movmtcd 
to the gimbal frame with a three point .suspension. The o*5tlcal bench is 
Separated from the balance of the LRS compox\ents by a thermal curtain to mlni- 
mlae radiative heat transfer. 

A typical ranging pass follows the following sequence^ The SIPS Is deployed 
and oriented for a prep.'iss cool-down period of currently unspecified dviration. 
After thernuil stabillKation has been achieved, the system is ready for the tiext 
ranging pass, and a wait mt>de is initiated. Approxinvately A5 minutes prior to 
the actual ranging pass, the LRS is selectively powered up and the SIPS is 
commanded to point to several stars to reduce the initial attitude reference 
uncertainty to acceptable limits. A programmed sequence of stellar sightings 
Is used to calibrate the ARS gyros. The Shuttle /Spacelab is then coninanded to 
the desired inertial attitude for the ranging pass, and the SIPS is oriented 
to the nominal inertial position for acquisition. 

Thirty seconds befoi*o the first target comes into view, the laser is energized 
and allowed to stabilize. Wlien the first target is expected to be in view, an 
acquisition sequence is initiated. A short target verify period follows initial 
target detection. The experiment then enters the normal ranging mode, which 
continues until the last target Is no longer in view. At this time, the I.RS is 
shut down, and the SIPS returns to the cool-down mode. 

Data accumulated during the pass is recorded for subsequent transmission to 
the ground. A typical Individual measurement record Includes both the time 
the pulse left and the time of arrival of the reflected pulse, gimbal angle data, 
a target label, and miscellaneous system parameters pertinent to data reduction, 

3.2,1 Functional Descrip tlon 

The major elements of the laser ranging system are shown in Figure 51. The 
major Interfaces are shown to illustrate the more significant function.s and 
data flow requirements. The clock subsystem issues a command sequence to the 
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laser Cransmltter subsystem to Initiate a pulse at, or near, an expected time, 
nominally at a 10 pps rate. The clock subsystem also determines the actual 
pulse departure time and the time of arrival of the reflected pulse. These 
data are forwarded to the data assembler and the navigation computer. The 
navigation computer either accepts or rejects the data point based on several 
test criteria. If accepted, the measurement is used to update the state 
vector via a Kalman filter technique. 

The navigation computer also generates a target range vector and range rate 
vector for the next pulse emission time and forwards this data to the ARS 
computer. The magnitude of the range vector is used to set the range gate In 
the ranging receiver. The range and range rate vectors are processed in the ARS 
computer to determine the current pointing error and the command sequence 
necessary to drive the SIPS gimbals to minimize the pointing error at the next 
pulse emission time. Hiis control data must be output at a relatively high 
rate to ensure smooth response of the SIPS glmbal control system. Note that 
the actual glmbal angles are needed by the navigation computer (and the data 
assembler) in order to account for the range bias due to offset from the Shuttle 
center of gravity. 

Also shown is a television camera and control to enhance the ability of the 
experiment controller to monitor and control the experiment operation. 

The laser ranging experiment is conducted during a number of discrete ranging 
opportunities, which occur when the ground track of the Shuttle passes suffi- 
ciently close to a target grid. A typical ranging pass begins about 45 minutes 
prior to the point of closest approach. The first step is to power up the 
clock, the attitude reference system, and the navigation computer. The clock 
is Initialized to the Shuttle GMT. The ARS Is initialized to the Shuttle 
reference, rotated through the SIPS gimbal angles. The SIPS is then commanded 
to a preprograimned series of stellar sightings to refine the ARS attitude 
reference and to calibrate the gyro drift and scale factor. At a suitable 
time prior to commencement of the ranging pass, the navigation computer deter- 
mines the desired Shuttle attitude, in stellar inertial coordinates, and 
commands the Shuttle to maneuver to and maintain this attitude. The SIPS is 
then commanded to the acquisition orientation, and the normal pointing mode 
is initiated. Thirty seconds prior to the predicted time for the first target 
to come into view, the laser is energized and allowed to thermally stabilize. 
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Th« ayatem then enters the Acquisition mode, which requires n sequential search 
about the nominal pointlnR direction in a preprograiimved search sequence. n\e 
acquisition strategy selected for this experiment Is a varlabK'^ dwell approach. 
At each of the search locations, the received signal Is compared to dual 
threshold criteria. The search would continue at each location until a posi- 
tive no target or a positive target In view decision would be nuide. This 
strategy Is maximally effective for the relatively slow moving SIPS platform. 

Once the target Is detected, a short dither sequence is performed to estlnutte 
the true target angular location and the system enters the ranging mode. 

In the ranging mode, the SIPS Is commanded to point toward each of the targets 
In the grid in a preprogriwttwed sequence. At each discrete target, ranging 
measurements are made for n successive pulses, where n Is a variable dependant 
on the size of the grid and the angular agility of the SIPS. Kach measurement 
group is tested to ensure validity, and used to update the state vector using 
a Kalman filter technique. This sequence continues until the last target in 
the grid Is no longer In view. A typical ranging pass lasts for about 260 
seconds . 

Tl»o post-pass mode is provided to ensure an orderly power down sequence and 
completion of all buffered fvmcclous. The post pass mode la completed when the 
last of the LRS components is powered dowit and the SIPS returns to the cold 
soak mode. 

3,2,2 Opto-Mecha«lcal Subsystem 

The OMS Is virtually identical to Che subsystem described for the pallet 
mounted experiment. Tl»e iiuijor difference is that it is not necessary to fold 
the transmit beam to the receiver centerline, since the glmballed mirror is 
not required* The laser output is passed through a frequency doubler and thou 
expanded to approxinately 2 cm diameter* The frequency doubled output Is 
reflected by a dlchrolc mirror, and delivered to the oxitput lens, which 
adjusts the beam to a 0,5 mrad divergence. The beam also passes through a 
partially reflective mirror which extracts approxlnuitoly one percent of the 
pulse energy. This portion of the pulse energy Is passed through a neutral 
density filter and reflected to the ranging detector for a departure time 
measurement. 
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Th© receive optics psth includes a bifurcatluR ralin'or, which allows a star 
sausor to be used to calibrate the receiver boresijjht* The central portion ot 
the receiver field of view passes thvoush the blCmcating irdrvor, and encounters 
two torque motor beam steering components, a narrow band filter, and a variable 
neutral density filter. Another bifurcating mirror limits the high speed 
detector field of view to '\*0.5 lurad, and reflects the balance of the field of 
view to a quadrant photomultiplier. 

Various alternatives for on-orbit calibration of transmitter and receiver 
alignment were considered. The most appropriate technique was Judged to be 
performed tix conjunction with a cooperative ground station. The ground station 
would illumluate the Shuttle. The QVMT output is used to detormlue the precise 
calibration of receiver boresight to ARS baseplate. The l.RS transmitter output 
would be detected by a cluster of calibrated receivers on the ground, and a 
boresight error determined. Tit Is data would then be transmitted to the l.RS via 
normal Shuttle comnumi cat ions, and used to bias the computed potutlug angles. 

A televi.slou camera could be added to enhance the capability of the experiment 
operator to monitor and contixil the experiment. The collecting power and resolu- 
tion of the receiver telescope would allow both the star sen.sor and the tele- 
vision c.araera to share iSO-30 division) the ovitev portion of the receiver field 
of view, .\lternatively , a torque motor be.am steering element could be employed 
to select either the star sensor or the televl.sion camera. 

Interfaces 

The experiment ii^terfaces are with the SIPS, and nltimtjtely with the Shuttle/ 
Spacalah. The hRS optical bench Is attached to the catiisttsr glmbal frame and 
thermally shielded as much as possible. Those components not nK>ni\ted on the 
optical bench are attached to a support structure behind the optical bench. 

The LRS data interface with the SIPS and the Shuttle/Spacel.'ih is through a 
standard RAH. .As presently planned, data will be transferred to the Shuttle 
for storage and subsequent transmission to the grovind. In the same manner as 
for the pallet mounted experinmnt, and will not be prtisex\ted again In this sec- 
tion. The total data load through the KAU will, however, be somewhat greater 
for this experiment, due to the need to transmit pointing coimmtnds to the SIPS 
computer anC to receive glmbal angle data necessary to eliminate the time 
varying range bias. 


:u3i 
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The electrical power budget shown In Section 3.1,4 Is the same except for dele- 
tion of the thermal control subsystem (provided by the SIPS). The energy 
budget, however, Is considerably reduced since the clock will not be operated 
continuously. Table 11 summarizes tlie energy consumption for this experi- 
ment. The prime power requirement for the power conditioner has been adjusted 
to account for deletion of the active thermal control subsystem. 
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TABLE 11 

POWER PROFILE ESTIMATE 



POWER VS 

ALIGN 
(45 MIN.) 

MODE 

PASS 

(280 SEC) 

POSI-PASS 
(5 MIN) 

ENERGY 
CONSUMED 
(10 PASSES) 
(KW-HRS) 

LASER PWR SUPPLY 

0 

39GW 

0 

0.3 

ARS SENSORS 

66W 

66W 

66W 

0.6 

POWER CONDITIONER 

13 7W 

13 7W 

137W 

1.25 

TOTAL 

2G3W 

593W 

203W 

2.15 


SYNTHESIS AND ANALYSIS OF 

PRECISE SPACEBORNE LASER RAN61NS SYSTEMS 


REPORT MDC El 729 

VOLUME II 

FEBURARY 1978 


A. Conclusions 

A laser ranging experiment has been defined and analyzed • and two laser 
ranging systems have been configured to accomplish the experiment. The experi'- 
ment has been shown to be capable of meeting the target grid reconstruction 
accuracy goal. A number of Important factors have been identified and docu- 
mented, such as the Importance of widely spaced out-rigger targets. The 
prototype system used to validate the experiment design was a one pps laser 
ranging system with a 10 cm rms accuracy, with relatively modest slew rate 
requirements. The two laser ranging systems employ a 10 pps laser, with rms 
accuracies In the 2 to 10 cm range, depending on environmental considerations. 
The pallet mounted experiment has considerably greater angular agility than 
the SIPS mounted experiment, since only the glmballed folding mirror must be 
moved. Both systems are Judged capable of meeting the assumed experiment 
objectives of determining the relative positions of five targets with respect 
to each other with an overall three dimensional reconstruction accuracy better 
than 2 cm rms; the pallet mounted experiment is capable of accomplishing more 
ambitious experiments. 
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APPENDIX A 


Coupling of geopotentlal terms with spacecraft motion. 

Let X be spacecraft position 
V be spacecraft velocity 
and h be geopotentlal terms. 


Then the equations of motion are 

X - V 

■ 

V - a(x.h) 
h - o 

where a represents the total acceleration acting on the spacecraft. Then 



and the state transition matrix is 


At 

J F(T)dT 


At/o l\ 

4 -) 


(: :) 


dT 


Here 4> - e ^ is the usual orbit state transition matrix. To 

obtain note that 4) and <Jj satisfy 

, ei>(4‘ 4 )(::>•»(: 0 

‘ \ 0 0 0 / 


This gives 


/ 


\ 
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-(t 


giving 




;/;( 0 ) » 0 


/ (ft) 


dT 


We assume to be constant (l.e* the geopotentials do not change appreciably 
over the state transition time Interval). If we cannot do this then numerical 
integration would have to be resorted to and this we want to avoid. This 
assumption allows 3a/ 3h to be factored from the Integrt;.^! 


^(t) - <i»(t) 


f <t*~^(i)dT/ 3a \ 

; \ ; 


■/ 


e- 


t 

■/ 

t 

■/ 


<l>(t) (Ji“^(T)dTl 3a 


<}>(t)(()(“T)dT 





80 that coupling terms can be obtained by simple quadratures from the space- 
craft state transition. 


SYNTHESIS AND ANALYSIS OF 

PRECISE SPACEBORNE LASER RANGING SYSTEMS 


REPORT MDC El 729 
FEBURARY 1978 


VOLUME II 


APPENDIX B 
POINTING COVARIANCE 


Following is the derivation of an expression for the variance of the pointing 
uncertainty for directing an Instrunwnt on the spacecraft with position vector 
r towards a reflector at p. Let S be the vector from F to p then 

P » P - r (1) 

I<ec e^ be a unit vector in the direction of P and denote vector magnitude by 
omission of overbars. We then have 


p - pe^ (2) 

It is desired to express the error in angular orientation of p due to the errors 
in r and p. The covariances of F and p are obtainable from the total state 
covariance matrix. The angular variation is given by a variation In e , From 


But 


5 — a 
P 


6 e + p6e 
P P P 


2 - - 
p - p • p 


giving 



(3) 

(A) 

(5) 


Inserting (5) Into (3) 

+ P<Se 

\P I p ^ 

or using (2) 


p6e^ - 5p - 


Taking the dot produced of (7) with Itself 

-2(Ip * 6p)^ 

® 6p ‘ Sp - Cfip • 6p)(5p 


ep) 


B-1 


( 6 ) 

(7) 

( 8 ) 




- I '"3 5 
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The angular variance of Interest Is 0^ 
(8) we obtain 


6p^) so taking expectations of 


E(6p - 6r) • (6p - 6r) 


“ • E((fip - 6r)(6p - Sr)) • (9) 

where the central factors of the last tern is a dyadic. We let the 3x3 
matrix A represent this dyadic then 

A-S -S -S +S (10) 

pp pr rp rr 

where S is the covariance of the reflector position components, S is the 
pp rr 

Covariance of spacecraft position and is the correlation between 

spacecraft and reflector position (induced by the estimation procedure). 
Expression (9) becomes 

(tr A - e*^Ae) (11) 


where tr A denotes the trace of A (l.e., trA - + ^22 ^ ^33^ e is a 3 x 1 

column vector representing to vector e^ (l)e., e contains the direction cosines 
of F)* Expression (11) is implemented in a subroutine. 
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APPENDIX C 

Suppose that one observation Is made on reflector //I from S and subsequently 

one observation on reflector 92 * Assume that the geometry of the observer 

reflector configuration does not change between observations and that the 

2 

observations are of range with a random error having variance Q » o and a 

2 ® 

systematic error having a variance Pj, * • W^e wish to find the 3x3 covari- 

ance matrix P of the differential position of reflector #1 with respect to 
reflector #2. 

Let the state vector be 

( 3 coordinates of reflector #1 ^ 

3 coordinates of reflector #2 

3 coordinates of observer 

Systematic error In measurement^/ 


10 X 1 

X • 


3x1 
rl 
3x1 
r2 

3x1 
rS 
1 

rb 


\ 


hJ 


3'X 3 

Let both reflectors have an aprlor position error covariance P and the 

3x3 ^ 

observer position error p , then If these are uncorrelated the apriori 


covariance^ of the 

state 

Is 



h 

0 

0 0 \ 

10' X 10 

r 


0 0 

P 

o 

r 

0 

-S 0 , 


\» 

0 

0 


( 1 ) 


1x1 


For a measurement of range from observer to reflector #1 we have y^^ 
determined by the expression 


(Yj - b)‘ 


(r^ - S)^(rj - S) 


where the superscript T denotes transpose of a matrix. Taking variations in 
this we obtain 


2(Yj - b)(6y^ - 6b) 


2(r^ - S)*(6,r^ - 6S ) 


giving 




(r^ - S) 


(6r^ - 6S) + 6b 
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Let a unit vector In the direction of reflector ifl from the observer be denoted 
by and similarly 62 for reflector #2. Then 



and 


6Y^ " - ^S) + fib 

1x3 1x3 1x3 1x1 

- 0 . 1 ) 

- 0, -ej) fix 

Denote 



" a 

- 




1 ) 


3X 


( 2 ) 


This will be needed In to determine the stato covariance after the observation. 
The Kalman update formula is 


10x10 10x10 10x10 10x1 1x10 10x10 10x1 1x1 

P, - P - P (M, P 


1x10 10x10 


1 ■„ fo 


<3) 


where P^ is the state covariance after incorporation of the observation on 
reflector ^ 1 . 

1x1 


Define + Q 


(A) 


Now from (1) and (2) 


P M,T 
o 1 


0 

0 

\o 


0 

0 

0 


0 

0 

"s 


■Alt 

k! 


\1, 




T‘^l\ 


"^s®li 




(5) 
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and 






* *l '■s'l 




+ Q 


and from (3) (factoring out the scalar ^ using (3) and noting that 
Pj^ are symmetric 




/**T 

0 

0 

\o 


T 

0 P^ 
0 0 


0 
0 




l/ZlP^e^e^ P^ 


For the observation of reflector #2 we have 




1. 


■ "l 

i-vJ 


\'b / 

T ® 

l/Zj! 


0 

0 

**s- 

0 

1/ZjP, 

we have 



0. -e^\, P^ 


-l/ZlP^eiP 


'1 S 1 1 -S l/2iPse,P, 


IxlO T I 

M- - (0. e_\ - eJ ,1) 


and 


with 


10x10 

P„ 




1. 

z. 


h”2\h‘ 


1x1 

^2 •Vi«2^Q 
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( 6 ) 



( 8 ) 


(9) 


( 10 ) 
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Procedlng as before 


P M 
12 




T 2 


\ 


~^S®2 ^s®2 ^b^ 

"b ^ i; ^ "b> 


®2 ^ T ®2 '*■ ®2 ^ S ®2 “ ®2 ^ S ®2 ^ b ^ 

^b “ ^b ^S®2 + P^) + Q 


We are not directly Interested in the covariance P 2 of the entire state vector 
but rather in the covariance of the separation p of the reflectors. Let 


3x1 3x1 3x1 

P - t2 - Ti 


then 


l£ 


3x10 3x3 3x3 3x3 3x1 

gjj « H = (-1, I, 0. 0) 

where I ■ / 1 0 o\ 

0 i 0 

0 0 1 

. / 

The covariance of the error in p after both observations Is 


( 12 ) 


3x3 

’■p 


3x10 10x10 10x3 


H 


H' 


HP^h'^ - 1/22 HJ’i”2^“2*’i“^ 


by referring to (9) , to evaluate this expression ve will do It piece by piece as 
follows . 
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From (7) and (12) 


HP^H* - 2P^ 




and from (7), (12) and (8) 

“A'' - V2 + * V 

Inserting these into (13) we obtain 

^ Vi"i% 

■ ■*' ^S®2 z!j^ ^®2 ^S®1 

1 T 

"p ■ ®T - •'i 


Let 1x1 A 


" ®1 ^S®2 ■*' 


(15) 


then 


P^ ” 2P^ “ ^T®1®1 


^P^(e2+|e^)02^ + ^e/)PT 


(16) 


Iteftae - Z^ - e/(Pj + PjXe^ + P^ + Q 

(17) 

^20 - *2^«I * ^>'2 

(18) 

then 

Z 2 


^2 • ^20 * ^ 

(19) 

by referring to (11). 
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A factor appearing in the middle of the last two terms of formula (16) Is 


11 * Zj '°2 Zj^ ® l ' '*2 * Zj ®1 ^ 


T T 

(e, +^ej(e, + ^ e/ 


2 Z^g“l 


"20 2 


10 

\ 


1 ^ "S 1 

^10 ^10 , 2 


"20 z^^y 


eie^ + 


1 . S 1 . T. T, 

e„e„ + s +eie-i ) 


2^2-2 ^10 
^20 - 2^0 ^20 - 




"20 


T 

e,e, + 


"10 


e«e„ + 


T T 

(e,e, + e.e, ) 


2 ‘ 2 ^^2^2 *" 2 '21 12 

W20-"S WZO-^S VZO-^S 


SO that the formulas (16), (15), (17) and (18) have the form which, as expected, is 
symmetric with respect to the two measurements 


F *> 2 P - P 

T T 


^20 V/ + SoV2^ + . 

2 ^ 

^10^20 " 


’ ®1 ^S®2 


Zio » e, (P^ + Pg)e^ + + Q 


Z20 - (P^ + Pg)e2 + Pj, + Q 


Special case: Both reflectors with same direction from the observer. 


In this case e^^ « 62 * e 


From (20) 


^10 “ ^20 “ ^0 


( 20 ) 


and 
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P « 2P - P ^^0 T 

ee\ 


2P„ -2P^ P„ 

T T - Z„ T 


with Z^ - Zg = (e’^(P^ + Pg)e + + Q) - (e’^’Pge + P^) 

“ e^^P^e + Q 

The minimum eigenvalue of P^ must obviously correspond to the direction of e and 
is given In this case by 

T„ 

T ® 

' V - ' ^ 

e P^e + Q 
T 

2 ® 

since Q » a and -= — — — < 1 

” + Q 

We have that the standard deviation of the psoition difference In the direction 
of the measurements Is 


< \/T o 

m 

as suspected. 

Suppose the reflector directions are separated by an angle 26. Set up a 
coordinate system with x axis bisecting the directions to the targets with the 
y axis in the plane of the observer and both reflectors. The unit vectors e^ 
and e^ are 

( cos 6 \ / cos 6 ' V 

sin I , ®2 * ( ^ I 

Denote the unit vector in the direction of the x-axis by e 


0 


Assume for simplicity, that In this coordinate system that all of the aprlor 
position covariances are diagonal. Thus, 
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TX 


Ty 


0 \ 
0 


L ^ 


P - 


sx 


0 

\0 


sy 


0 \ 
0 

■’sz/ 


Q « er 
m 

T. 2 


Referring to (20), 


2 2 2 2 2 
Z„ “ CT„ eos 5 - CT„ sin 6 + ct 
S Sx Sy b 

2 

^0 “ ^10 ■ ^20 " ■*“ "^Sx^^ ^ + ^^y )aln^ « + 

From symmetry we know that the minimum eigenvalue of P corresponds to the 

T 

direction e so we wish to evaluate e 

T 2 T T 2 

Now e P^e « 0 ^^ , e P^e^^ - e^ P^e » 0 ^ cos & 


and 


2 ^20°Tx* ^ * ho'^Tx 

0 ' ^ n — 

^10 20 


2 4 2 _ 

- “tx ““ ' z/ - 2," 


^^Tx „ 2 2 ,, 

Zq - Zg ^^0 ~ ~ Tx ^ 


20 , 


, 2 , - 2 . ^ 2 ^ 2 

2 <V ^°Sy ^ *’ * °n 

Tx 2 2 2 2 2 2 

0_ cos 5 + (0„ + 20 _ )sln 5+0 

Tx Ty ■ Sy m 


This Is the desired formulav It shows how the separation angle 6 

2 

allows the observer position variance ) (normal to the mean line of sight 

and in the plane of the observer and targets) to degrade the experiment since 

2 2 2 
0 g sin 5 can easily be of the same order of magnitude as 0 
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